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ABSTRACT 


The  equivalent  circuit  representation  of  the  WR-90  waveguide  resonant 
circuit  with  post  mounting  of  devices  is  used  to  analyze  the  HF  loads 
presented  to  X-band  Guim-effect  devices.  Experimental  data  for  a  pulsed 
oscillator  is  analyzed  over  the  frequency  range  of  8  to  12  GHz.  Domain  mode 
behavior  is  found. 

Calculated  bandwidths  of  the  order  of  an  octave  are  shown  for 
stabilized  transferred-electron  amplifiers  using  short  crystal,  lengths. 

Effects  of  lattice  temperature  on  the  behavior  of  these  devices  are  described. 
The  possibility  of  improving  the  dc  to  RF  conversion  efficiency  by  a  lattice 
temperature  gradient,  or  by  a  doping  gradient  along  the  length  of  the  GaAs 
crystal  is  discussed. 

The  frequency  modulation  sensitivities  of  a  Gunn-effect  device 
mounted  in  a  coaxial  circuit  is  measured  for  different  lattice  temperatures 
and  bias  values.  The  rapid  changes  in  modulation  sensitivity  are  shown 
to  be  due  mainly  to  device  changes  and  not  RF  circuit  properties,  m 
addition,  the  circuit  parameters  that  limit  the  electronic  tuning  rate  of 
the  device  in  a  waveguide  circuit  are  found. 

Equivalent  circuits  for  the  diode  package  transformation  in  50 -SI 
microstrip  line  and  in  ridged  waveguide  are  presented,  the  method  for  their 
determination  is  discussed  and  admittance  characteristics  for  GaAs  and  Si 
diodes  biased  at  1000  A/cm2  are  given. 

Further  development  of  the  computer  program  modeling  of  the  IMPATT 
mode  of  the  a^ilanche-diode  oscillator  has  been  carried  out  which  has 
resulted  in  a  significant  improvement  of  its  rate  of  convergence.  A  program 
to  calculate  independently  the  small-signal  admittance  for  an  arbitrary 
doping  profile  has  also  been  written.  A  description  is  given  of  the 
approach  to  be  taxen  to  permit  modeling  of  the  TRAPATT  mode. 

Some  representative  numerical  results  obtained  from  a  large-signal 
IMPATT  diode  computer  analysis  are  presented  and  discussed. 


An  improved  diode  fabrication  process  is  described.  A  group  of 
diodes  with  low  thermal  resistances has  been  obtained.  IMPATT  diodes 
fabricated  with  the  process  have  been  tested  with  powers  in  the  100  to 
600  mW  range  (CW)  and  with  frequencies  between  6.5  and  8.5  GHz. 

A  variable  step  size  numerical  integration  scheme  has  been  implemented 
for  calculating  the  collision  integrals  in  the  dc  analysis  of  germanium. 
Modifications  are  being  made  in  the  computer  program  so  that  this  scheme 
does  not  take  excessive  computer  time. 


TABLE  OF  CONTENTS 


Page 


1.  GENERAL  INTRODUCTION  1 

2.  HIGH-POWER  GUNN-EFFECT  DEVICES  1 

2.1  Introduction  1 

2.2  RF  Circuit  Characteristics  1 

2.3  Experiments  with  a  Pulsed  TE  Device  in  tile  Waveguide  Circuit  8 

2.4  Conclusions  17 

2.5  Program  for  the  Next  Quarter  l8 

3*  GUNN-EFFECT  AMPLIFIERS  l8 

3 . 1  Intro  duct ion  l8 

3*2  Bandwidth  Considerations  19 

30  Field  Dependence  of  Diffusion  Coefficient  and  Its  Importance 

in  Determining  the  do  Characteristics  Accurately  22 

3.4  Effects  of  Lattice  Temperature  22 

30  Variation  of  the  Lattice  Temperature  Along  the  Length  of 

the  GaAs  Device  30 

3.6  Variation  in  Donor  Density  Along  the  Length  of  a  GaAs 

Crystal  32 

3*7  Conclusions  33 

3.8  Program  for  the  Next  Quarter  53 


-v- 


Page 


4.  MOIXJLATION  PROPERTIES  OF  GUNN-EFFECT  DEVICES  34 

4.1  Introduction  34 

4.2  Experimental  Investigation  of  Modulation  Properties  of  a 

Gunn-Effect  Device  at  Different  Temperatures  34 

4.3  Device  Operation  in  Waveguide  and  Coaxial  Circuits  45 

4.4  Waveguide  Mounting  of  Devices  49 

4.5  Conclusions  50 

4.6  Program  for  the  Next  Quarter  51 

5.  AVALANCHE-DIODE  AMPLIFIERS  51 

5.1  Introduction  51 

5.2  Determination  of  Immittance  Transformations  Characterizing  a 

Diode  Package  and  Mount  52 

5.3  Package  Transformation  in  a  Ridged  Waveguide  57 

5.4  Admittance  Characteristics  of  Si  and  GaAs  Diodes  57 

5.5  Program  for  the  Next  Quarter  60 

6.  HIGH-EFFICIENCY  AVALANCHE  DIODES  60 

6.1  Introduction  60 

6.2  IMPATT  Mode  Computer  Program  63 

6.3  TRAPATT  Mode  Program  66 

6.4  Conclusions  67 

6.5  Program  for  the  Next  Quarter  6 7 

7.  NONLINEAR  OPERATING  CHARACTERISTICS  OF  IMPATT  DIODES  67 

7.1  Introduction  67 

7.2  Nonpunch -Through  Silicon  Diode  67 


’•V«" 


Page 

7*3  Silicon  Abrupt -Junction  Diodes  with  5  x  10lS  Background 

Doping  75 

7-3 •!  Small -Signal  Analysis  76 

7.3*2  Large-Signal  Analysis  8l 

7.4  Millimeter-Wave  Silicon  Diodes  90 

7.5  Conclusions  and  Program  for  the  Next  Quarter  98 

8.  SOLID-STATE  DEVICE  FABRICATION  98 

8.1  Introduction  98 

8.2  Diode  Fabrication  99 

8.3  Experimental  Results  '  ',1 

8.4  Conclusions  106 

8.5  Program  for  the  Next  Quarter  106 

9.  INSTABILITIES  IN  GERMANIUM  106 

9*1  Introduction  106 

9.2  Modification.'  in  the  Simpson's  Integration  Procedure  107 

9.3  Results  108 

9*4  Program  for  the  Next  Quarter 


109 


LIST  OF  ILLUSTRATIONS 


Figure 


2.1  Load  Resistance  Presented  to  the  TE  Device  in  the  Waveguide 
Circuit  as  a  Function  of  Frequency.  Short  Plane  Is  2.354  cm 

from  the  Post.  3 

2.2  Load  Reactance  Presented  to  the  TE  Device  as  a  Function  of 
Frequency  in  the  Waveguide  Circuit.  Short  Plane  Is  2.35^  cm 

from  the  Post.  4 

2.3  Load  Admittance  Presented  to  the  TE  Device  in  the  Waveguide 
Circuit  as  a  Function  of  Frequency.  Short  Plane  Is  2.354  cm 

from  the  Post.  5 

2.4  Lo  .d  Admittance  Presented  to  a  TE  Device  in  the  Waveguide 
Circuit  as  a  Fune'-.ion  of  Frequency.  Short  Plane  Is  1.755  cm 

frc  i  the  Post.  7 

2.5  Load  Conductance  as  a  Function  of  Resonant  Frequency  for 

Matched  Waveguide  Circuit  for  Two  Values  of  Load  Susceptance.  9 

2.6  Pulsed  RF  Power  Output  as  a  Function  of  Frequency  for  TE 

Device  M61-24  in  the  Waveguide  Circuit.  (Vg  =  15  V)  10 

2.7  RF  Voltage,  Load  Conductance  and  Susceptive  as  a  Function  of 

Frequency  of  Oscillation  for  Pulsed  TE  Device  M61-24  in  the 
Waveguide  Circuit  Without  a  Tuner  ind  with  15  V  Bias.  11 

2.3  Load  Conductance  and  Susceptance  vs.  Frequency  for  TE  Device 
M61-24  for  the  Matched  Load  Case  and  for  the  Best  RF  Power 
Condition  Using  an  RF  Tuner  in  Load  Lino.  (Vg  =  15  V)  13 

2.9  Admittance,  Capacitance  and  RF  Voltage  as  a  Function  of  RF 

Frequency  for  TE  Device  M61-24  in  the  Waveguide  Circuit  with 
a  Tuner  for  Best  Output  Power.  (Vg  =  15  V)  l4 


Figure  Page 

2.10  Admittance,  Capacitance  and  EF  Voltage  as  a  Function  of  RF 
Frequency  for  TE  Device  M61-24  in  the  Waveguide  Circuit  with 

a  Tuner  for  Best  Output  Power.  (Vg  =  20  V)  15 

2.11  RF  Power  Output,  Frequency  of  Oscillation,  EF  Voltage,  Load 
Conductance  and  Susceptance  as  a  Function  of  Bias  Voltage 
for  TE  Device  M6l-2^  in  the  Waveguide  Circuit  Without  a 


Tuner  and  for  a  Post-to-Short  Distance  of  2.35^  cm.  16 

3.1  Impedance  as  a  Function  of  Frequency  at  300°K  with  Diffusion 
Effects  Included.  (nQ  =  5  x  10l4/cm3,  Length  =  l4  um, 

J  =  1.1  and  Area  =  9*8  x  10-4  cm2)  20 

o 

3.2  Impedance  as  a  Function  of  Frequency  at  300°K  with  Diffusion 
Effects  Included.  (nQ  =  5  x  10l4/cm3,  Length  =  10  pm, 

J  -  1.1  and  Area  =  7  x  10  4  cm2)  21 

o 

3 '3  Average  Carrier  Density  and  Diffusion  Coefficient  vs.  Distance 
from  the  Cathode  at  300°K  with  Diffusion  Effects  Included. 

(n  =  5  x  1014/cm3,  J  =  1.1,  J  =  0.001  and  Frequency  = 
o  °  11 

10  GHz)  23 

3-^  Velocity -Field  Relationships  Used  for  Small-Signal  Diffusion- 

Free  Study.  25 


3.5  Diffusion-Free,  Small-Signal  Calculations  of  Maximum  Negative 
Resistance  and  Length  of  Active  (i^)  and  Passive  (L&)  Regions 

in  a  Device  10  urn  Long  with  nQ  =  5  x  10l4/cm3  and  Area  = 

7  x  10  4  cm2.  Velocity-Field  Relations  of  Fig-  3*^  Were  Used.  26 

3«6  Static  Electric  Field  vs.  Distance  from  the  Cathode  for  a 

Small-Signal  Case  at  Different  Lattice  Temperatures  with 

Diffusion  Effects  Included.  (J  -  1.1,  'j  =  0.001. 

o  11 

Frequency  =  10  GHz  and  nQ  =  5  x  10l4/cm3)  28 

3.7  Impedance  vs.  Frequency  at  Different  Lattice  Temperatures  with 
Diffusion  Effects  Included.  (no  =  5  x  10l4/cm3,  Bias  Voltage  = 

lb  V,  Length  =  10.1  (am,  Area  =  7  x  lO-4  cm2  and  J  =  0.I5)  29 

11 

3.8  Dc  Device  Current  vs.  Bias  Voltage  at  Two  Different  Lattice 
Temperatures.  (Length  =  10  (am.  Area  =  3.5  x  10~4  cm2  and 

nQ  =  5  x  10l4/cm3)  31 


-ix- 


Figure  pa go 

4.1  Experimental  Setup.  35 

4.2  I-V  Charrcteristics  of  Device  KB-10B-105  for  Heat  Sink 

Temperature  at  8.5°C.  The  Low-Field  Resistance  Is  8  Q.  36 

4.3  I-V  Characteristics  of  Device  HB-10B-105  for  Heat  Sink 
Temperature  at  -38°C.  Enlarged  Low -Field  Characteristic 

and  Low-Field  Resistance  Are  Drawn  Under  the  Curve.  37 

4.4  I-V  Cijuracteristics  of  Device  HB-10B-105  for  Heat  Sink 
Temperature  at  -6l°C.  ESilargea  Low-Field  Characteristic 

and  Low-Field  Resistance  Are  Drawn  Under  the  Curve.  38 

4.5  Frequency  Modulation  Sensitivity  for  Heat  Sink  Temperatures 
at  -38°C  (a)  as  a  Function  of  Peak  Modulation  Voltage  and 

(b)  as  a  Function  of  Modulation  Frequency.  4o 

4.6  Frequency  Modulation  Sensitivity  vs.  Bias  Voltage  for 

Different  Heat  Sink  Temperatures.  4l 

4.7  Frequency  Modulation  Sensitivity  vs.  Device  Average 

Temperature.  42 

4.8  Amplitude  Modulation  Sensitivity  for  Heat  Sink  Temperature 
at  8.5°C  and  RF  Frequency  at  8.56  GHz.  (a)  As  a  Function 
of  Modulation  Frequency  when  the  Modulation  Voltage  Is 

at  0.2  V  and  (b)  as  a  Function  of  Modulation  Voltage.  43 

4.9  Amplitude  Modulation  Sensitivity  for  Heak  Sink  Temperature 

at  -6l°Cand  RF  Frequency  =  8.68  GHz.  (a)  For  Peak  Modulation 
Voltage  at  0.6  V  as  a  Function  of  Modulation  Frequency  and 
(b)  as  a  Function  of  Peak  Modulation  Voltage  at  Different 
Bias  Voltages.  44 

4.10  (a)  Calculated  Coaxial  Cavity  Circuit  and  Device  Lines  for 
the  Experiment  and  (b)  Enlarged  Figure  near  Point  2  in 

Fig.  4,10a.  46 

4.11  (a)  Typical  Waveguide  Cavity  Admittance  and  Negative  Device 
Admittance  for  a  Constant  Bias  Voltage  and  (b)  Enlarged 
Figure  near  Point  2  in  Fig.  4.11a  and  Device  Lines  for 

Different  Bias  Voltages.  47 

5.1  Measured  Impedance  for  a  Diode  in  an  023  Package  Terminating  a 

50-fi  Microstrip  Line.  Inset  Shows  Postulated  Pole-Zero 
Characterization.  55 

5.2  Assumed  Equivalent  Circuit.  55 


-x- 


Figure  „age 

5«3  Measured  and  Calculated  Impedances  for  a  Diode  in  an  023 

Package  Mounted  in  Ridged  Waveguide.  58 

5.4  Equivalent  Circuit  for  a  Type  023  Package  in  Ridged 

Waveguide.  59 

5.5  Comparison  of  Admittance  Characteristics  for  Unpackaged 

GaAs  and  Si  Diodes  Biased  at  1000  1./ cm2.  6l 

5.6  Terminal  Impedance  for  a  Packaged  (Type  023)  GaAs  Diode 
in  a  50-ft  Microstrip  Line  as  a  Function  of  HP  Voltage 

on  the  Wafer  and  Frequency.  (Jq  =  1000  A/cm2)  62 

6.1  Effect  of  Mesh  Size.  64 

7.1  Dependence  of  Drift  Velocity  upon  Electric  Field  for 

Electrons  and  Holes  in  Si.  69 

7.2  Field,  Particle  Current  and  Hole  Density  Profiles  at 
cot  =  90  Degrees  for  a  9-pm  Si  Diode.  (Operating  Point: 

Vw  =  1  V,  Frequency  =  12  GHz  and  Jdc  =  500  A/cm2)  71 

7.3  Waveforms  of  Injected  and  Induced  Current  Density  for  a  9-|_un 
Si  Diode.  (Operating  Point:  Vgp  =  1  V,  Frequency  =  12  GHz 

and  Jdc  =  500  A/cm2)  73 

7.4  Waveforms  of  Injected  and  Inaoced  Current  Density  for  a 

5-um  Si  Diode.  (Operating  Point:  V™  =  1  V,  Frequency  = 

Kfc 

12  GHz  and  Jdc  =  500  A/cm2)  74 

7.5  Small-Signal  Admittance  Determined  from  the  Exact  Analysis 
for  Complementary  Si  Diodes.  (Background  Doping  =  5  x  1015 

cm  3  and  Width  =  5  |im)  77 

'”.6  Comparison  of  Exact  Small-Signal  Adraitt|nc£  with  the  Approxi¬ 
mate  Read  Diode  Expression  for  a  5-um  p  nn  Si  Diode  Biased 
at  "00  A/cm2,  'avalanche  Resonant  Frequency  =7*1  GHz  and 
ij l  -  Fractional  Drift  Length)  79 


7.7  Device  Admittance  as  a  Function  of  RF  Level  for  the  Comple¬ 
mentary  Si  Diodes.  (Background  Doping  =  5  x  1015  cm  3  and 
Width  =  5  um.  Operating  Point:  Frequency  =  l4  GHz  and 
o'dc  =  1200  A/cm2) 


82 


«  ■*■*.  M.J  <■*  'I  ^ 


Figure 


7.8  Comparison  of  Induced  Current  Waveforms  Determined  by  the 
Large-Signal  Models  for  a  5 -Mm  n+pp+,  5  x  1015  cm"3  Diode. 

(Operating  Point:  V^j,  =  50  V,  Frequency  *  l4  GHz  and 

Jdc  =  1200  A/cm2)  84 

7.9  Particle  Current  Profiles  a4  Several  Times  in  the  EF  Cycle 
for  a  5 -Mm  n+pp+,  5  x  1015  cm  3  Diode,  (operating  Point: 

Vrf  =  ^  Frequency  =  l4  GHz  and  Jdc  =  1200  A/cm2)  85 

7.10  Comparison  of  Particle  Current  Profiles  at  ojt  =  180  Degrees 
Determined  by  the  Large-Signal  Models  for  a  5-Mm  n+pp+, 

5  x  1015  cm  3  Diode.  (Operating  Point:  V'  =  50  V, 

Kr 

Frequency  =  l4  GHz  and  J,  =  1200  A/cm2)  86 

dc 

7.11  Comparison  of  Induced  Current  Waveforms  Determined  by  the 
Large-Signal  Models  for  a  5-Mm  p+nn+,  5  x  1015  cm  3  Diode. 

(Operating  Point:  V™,  =  50  V,  Frequency  =  l4  GHz  and 

Jdc  =  1200  A/cm2)  87 

7.12  Dependence  of  large- Signal  Admittance  upon  Effective  Avalanche 
Width  for  a  Diode  with  '’’otal  Width  =  5  Mm.  (Operating  Point: 

Vgp  =  50  V,  Frequency  =  l4  GHz  and  Jdc  -  1200  A/cm2)  89 

7.15  Comparison  of  Exact  Small-Signal  Admittance  with  the  Approxi¬ 
mate  Read-Diode  Expression  for  a  0.5-Mm  p+nn+  Si  Diode  Biased 
at  104  A/cm2.  (Avalanche  Resonant  Frequency  =  58  GHz  and 
Relative  Avalanche  Width  =  50  Percent)  92 

7.14  Device  Admittance  as  a  Function  of  Frequency  and  RF  Voltage 

for  a  0.5-Mm  p+nn+  Si  Diode  Biased  at  104  A/cm2.  95 

7.15  Power  and  Efficiency  as  a  Function  of  RF  Level  for  a  0.5-Mm 
p4nn+  Si  Diode,  (operating  Point:  Frequency  =  3o  GHZ  and 

Jdc  =  104  A/cm2)  95 

7.16  Comparison  of  Induced  Current  Waveforms  at  Several  RF  Levels 
for  a  0.5-Mm  p+nn+  Si  Diode.  (Operating  Point:  Frequency  = 

80  GHz  and  Jdc  =  104  A/cm2)  90 

7.17  Electric  Field  and  Electron  Density  Profiles  at  Several  Times 
in  the  RF  Cycle  for  a  0.5-Mm  p+nn+  Si  Diode.  (Operating  Point: 

Vpp  =  l4  V,  Frequency  =  80  GHz  and  T.  -  104  A/cm2)  97 


-xii- 


Figure 

8.1  Power  Output  and  Efficiency  as  a  Function  of  Frequency  for 
Various  Bias  Current  Levels.  Diode  1  Fabricated  from 
Wafer  D. 

8.2  Power  Output  and  Efficiency  as  a  Function  of  Frequency  for 
Various  Bias  Current  Levels.  Diode  2  Fabricated  from 
Wafer  D. 


LIST  OF  TAELES 


Table  Page 

4.1  Gunn  Device  Behavior  for  Several  Bias  Voltages.  48 

5-1  Equivalent  Circuit  Element  Values.  56 

7.1  Large-Signal  Data  for  Complementary  Diodes.  90 

8.1  Diodes  from  Wafer  C  (Double  Epitaxial  Material).  102 


-xiv- 


PUBLICATIONS  DURING  THE  LAST  QUARTER 


W.  E.  Schroeder  and  G.  I.  Haddad,  "The  Effect  of  Temperature  on  the 
Operation  of  an  IMPATT  Diode,"  Proc.  I EEE  (Correspondence),  vol.  59, 

No.  8,  pp.  1242-1244,  August  1971. 

W.  E.  Schroeder  and  G.  I.  Haddad,  "Avalanche  Region  Width  in  Various 
Structures  of  IMPATT  Diodes,"  Proc.  IEEE  (Correspondence),  vol.  59, 

No.  6  pp.  1245-1248,  August  1971. 

W.  E.  Schroeder,  C.  M.  Lee,  R.  J.  Lomax  and  G.  I.  Haddad,  "Avalanche 
Region  Width  in  Various  Structures  of  IMPATT  Diodes, "  Presented  at 
the  1971  Device  Research  Conf.,  Ann  Arbor,  Mich.,  June  28— July  1,  1971. 

A.  K.  Talwar  and  W.  R.  Curtice,  "A  Study  of  Stabilized  Transferred- 
Electron  Amplifiers,"  Presented  at  the  1971  Cornell  Conf.  cn  High 
Frequency  Generation  and  Amplification:  De/ices  and  Applications, 
Ithaca,  New  York,  August  1971. 

R.  W.  Laton  and  G.  I.  Haddad,  "The  ’'.x'ects  of  Doping  Profile  on 
Reflection  -Type  IMPATT  Diode  Amplifiers, "  Presented  at  the  1971 
European  Microwave  Conf.,  Stockholm,  Sweden,  August  1971. 


-xv- 


MICROWAVE  SOLID-STATE  DEVICE  AND  CIRCUIT  STUDIES 


1.  General  introduction  (G.  I.  Haddad) 

«pu»  research  under  this  program  is  concerned  with  investigation  and 
utilization  of  various  phenomena  and  techniques  for  the  generation,  amplifi¬ 
cation  and  detection  of  electromagnetic  energy  at  microwave  frequencies.  The 
work  is  oriented  mainly  toward  an  investigation  of  the  basic  properties  of 
avalanche  transit-time  and  transferred-electron  devices  and  the  incorporation 
of  these  devices  in  various  types  of  circuits  including  microstrip. 

Several  tasks  were  active  during  this  period.  The  status  of  each  and 
the  work  planned  for  the  next  period  are  discussed  in  detail  in  the 
following  sections  of  this  report. 

2.  High-Power  Gunn-Effect  Devices  (W.  R.  Curtice) 

2.1  Introduction.  Many  previous  experiments  with  pulsed  and  CW 
Gunn-effect  (or  transferred-electron)  oscillators  produced  data  incon¬ 
sistent  with  domain  mode  theory.  The  most  puzzling  aspect  was  the  inability 
to  tune  devices  over  a  wide  frequency  range  with  good  efficiency  in  the 
Q-mode  (quenched-domain  mode).  Part  of  this  difficulty  was  due  to  HF 
circuit  problems,  that  is,  the  inability  to  present  proper  impedance  to 

the  TE  device.  The  waveguide  circuit  analysis  and  pulsed  oscillator 
characteristics  described  in  this  report  will  enable  a  better  understanding 
of  these  frequency  tuning  problems. 

2.2  RF  Circuit  Characteristics .  The  equivalent  circuit  represen¬ 
tation  for  the  WR-90  waveguide  circuit  with  post  mounting  of  devices  was 
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presented  in  the  last  quarterly  progress  report.  The  physical  circuit 
consists  of  a  0.1-inch  diameter  post  centered  in  the  waveguide  and 
touching  a  TE  device  at  one  wall.  An  adjustable  waveguide  short  is 
behind  the  post.  An  BF  bypass  capacitor  exists  at  one  end  of  the  post 
so  that  device  bias  can  be  applied  through  the  post.  The  package 
parameters  of  the  TE  device  are  included  in  the  equivalent  circuit. 

Figures  2.1  and  2.2  show  the  calculated  load  resistance  and 
reactance  presented  to  the  TE  device  for  a  shorting  plane  position  2.354  cm 
from  the  post  and  for  three  different  waveguide  load  conditions.  ZT  is 

li 

the  waveguide  load  at  the  plane  of  the  post.  The  load  conditions  are 
Zq  (matched),  l/2  Zq  and  2  Zq.  The  second  and  third  load  conditions  can 
be  obtained  with  a  slj de-screw  tuner  of  VSWR  equal  to  2.0  at  a  proper 
distance  from  the  post.  Figure  2.2  shows  that  for  a  matched  load  the 
circuit  exhibits  a  series  resonance  at  8.802  GHz  and  a  parallel  resonance 
at  9*062  GHz.  The  circuit  is  inductive  only  between  these  frequencies 
in  X-band.  Oscillation  with  the  capacitive  TE  device  M61-24  occurs 
around  9*01  GHz.  At  this  frequency  both  load  resistance  and  reactance 
change  considerably  with  change  of  the  external  real  load  value. 

Figure  2.3  shows  the  load  admittance  for  the  same  circuit  and 
shorting  plane  position  and  for  the  same  external  loads.  The  admittance 
presentation  has  been  found  to  be  more  useful  for  TE  device  analysis. 

Note  that  oscillation  occurs  closest  to  the  frequency  of  parallel  circuit 
resonance. 

Figure  2.3  can  be  used  to  evaluate  the  circuit  factor  R0/Q0  from 


the  equation 


FREQUENCY,  GHz 


FIG.  2.1  LOAD  RESISTANCE  PRESENTED  TO  THE  TE  DEVICE  IN  THE  WAVEGUIDE 
CIRCUIT  AS  A  FUNCTION  OF  FREQUENCY.  SHORT  PLANE  IS  2.^b  cm 
FROM  THE  POST. 


FIG.  2.2  LCAD  REACTANCE  PRESENTED  TO  THE  TE  DEVICE  AS  A  FUNCTION  OF 

FREQUENCY  IN  THE  WAVEGUIDE  CIRCUIT.  SHORT  PLANE  IS  2.354  cm 
FROM  THE  POST. 
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LOAD  CONDUCTANCE  AND  SUSCEPTANCE,  mmho 


8.6  8.7  88  8.9  9.0  9.1  92  9.3 

FREQUENCY,  GH* 


j  .  FIG.  2.3  LOAD  ADMITTANCE  PRESENTED  TO  THE  TE  DEVICE  IN  THE  WAVEGUIDE 

.  CIRCUIT  AS  A  FUNCTION  OF  FREQUENCY.  SHORT  PLANE  IS  2.35*+  cm 

I  . 

FROM  THE  POST. 
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The  Rq/Qo  factor  for  the  matched,  circuit  is  0.99^  fi  while  it  is  1.55  ft  and 
0.868  fi  for  loads  of  l/2  Zq  and  2  Zq,  respectively.  The  external  Q  factor 
Qx  was  evaluated  for  the  matched  load  condition  and  found  to  he  200.  The 
load  conductance  estimate  used  in  earlier  work  of  (R  /Q  •  Q  )  1  « 

0.005  mho  (at  9 »0I  GHz)  is  only  a  little  different  from  the  actual  load 
conductance  of  0.006  mho.  However,  the  load  conductance  changes  so  rapidly 
with  frequency  that  this  method  of  calculation  is  not  reliable.  An  ideal 
parallel  resonant  circuit  would  have  nearly  constant  (device)  load 
conductance  around  resonance. 

Oscillation  ceases  for  low  bias  on  the  TE  devices  and  the  reason 
for  this  is  the  real  loading  change.  As  the  bias  is  lowered,  the  device 
susceptance  increases  which,  as  shown  in  Fig.  2.3,  lowers  the  frequency  of 
oscillation.  The  real  load  increases  as  the  frequency  is  lowered  and 
finally  becomes  too  large  to  sustain  the  circuit  controlled  oscillation. 

Calculations  for  matched  load  conditions  and  for  a  closer  short 
plane  position  show  much  larger  load  conductance  values  for  the  same  values 
of  load  susceptance.  Thus  without  an  external  tuner,  the  device  real  load 
is  greatly  increased  for  higher  frequency  operation.  Figure  2.b  sliows 
admittance  calculations  for  a  closer  short  plane  position  and  matched  load. 
It  is  seen  that  load  conductance  is  increased  at  parallel  resonance  and 
the  separation  between  series  and  parallel  resonant  frequencies  has  also 
increased  greatly.  The  R 0/Qp  value  at  10.5  GHz  is  i.-»5  which  is  not 
greatly  different  from  the  value  for  lower  resonant  frequencies. 
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f 

I 

3  FIG.  2.4  LOAD  ADMITTANCE  PRESENTED  TO  A  TE  DEVICE  IN  THE  WAVEGUIDE  CIRCUIT 

\ . 

(  AS  A  FUNCTION  OF  FREQUENCY.  SHORT  PLANE  IS  1.735  cm  FROM  THE 

POST. 
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Figure  2.5  shows  how  the  real  load  increases  with  increase  of  the 


circuit  resonant  frequency  (by  decrease  of  the  post-to-short  plane 
separation).  It  is  clear  from  this  data  that  for  the  same  values  of 
circuit  susceptance,  the  load  conductance  increases  strongly  with  increase 
of  resonant  frequency.  This  is  the  reason  for  the  poor  TE  device  perfor¬ 
mance  without  a  tuner  at  high  frequencies.  The  increased  loading  causes  a 
power  decrease  with  frequency  and  does  not  allow  oscillation  above  11  GHz. 

2.3  Experiments  with  a  Pulsed  TE  Device  in  the  Waveguide  Circuit. 
Preliminary  tests  of  TE  device  M61-24  were  described  in  the  last  quarterly 
progress  report.  The  device  has  the  following  parameters: 

n  =  2  x  1015/cm3, 
o 

l  =  13-7  pm, 

=  0.75  m2/V-s, 

area  «  10  4  cm2, 

tth  "  kA5v- 

I, =  0.525  A  (low  duty  operation)  and 

In 

Rlcw  field  =  5,52 

The  pulsed  power  output  obtained  'or  TE  device  M61-24  in  the 
waveguide  circuit  is  shown  in  Fig.  2.6  for  operation  with  a  slide-screw 
tuner  and  for  a  matched  load.  For  the  matched  load  condition,  the  device 
admittance  can  re  calculated  from  the  frequency  of  oscillation  using  the 
known  equivalent  circuit.  The  RF  voltage  can  be  obtained  from  the  RF  power 
value  and  the  real  load  presented  to  the  device.  These  results  are  presented 
in  Fig.  2.7  for  operation  at  15  V  bias.  The  increased  device  lo  ding  with 


frequency  not  only  reduces  the  RF  power  but  also  reduces  the  RF  voltage. 


o 


PIG.  2.5  LOAD  Cr^jJUCTANCE  AS  A  FUNCTION  OF  RESONANT  FREQUENCY  FOR  MATCHED 
WAVEGUIDE  CIRCUIT  FOR  TWO  VALUES  OF  LOAD  SU^CEPTANCE. 
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(b)  SLICE-SCREW  TUNER  USED  TO  MAXIMIZE  OUTPUT  POWER 

FIG.  2.6  PULSED  RF  POWER  OUTPUT  AS  A  FUNCTION  OF  FREQUENCY  FOR  TE  DEVICE 

M61-24  IN  the  WAVEGUIDE  CIRCUIT.  (V.  =  15  V) 

B 
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LOAD  CONDUCTANCE  AND  PEAK  RF 

NEGATIVE  SUSCEPTANCE,  mmho  VOLTAGE, 


CONDUCTANCE 


Operation  above  9*25  GHz  occurs  with  KF  voltages  insufficient  for 
single-frequency  Q-mode.  In  fact,  any  EF  voltage  value  between  2  and  13  V 
seems  to  be  allowable. 

The  load  presented  to  the  TE  device  with  and  without  a  tuner  are 
shown  in  Fig.  2.8  for  15  V  bias.  The  matched  load  results  are  slightly 
different  from  the  values  in  Fig.  2.7  because  a  correction  was  made  for  a 
slight  mismatch  (VSWR  <1.2)  in  the  load  line  which  was  found  to  be  present 
during  the  tests.  For  measurements  with  a  tuner,  the  actual  waveguide 
load  admittance  was  measured  after  adjustment  of  the  tuner  for  best  EF 
power  output. 

Notice  that  the  load  susceptance  is  nearly  the  same  for  the  two 
cases  in  Fig.  2.8,  but  that  the  load  conductance  is  significantly  lower 
at  high  frequencies  for  operation  with  the  tuner.  The  EF  voltage  differences 
existing  in  the  two  cases  mainly  affect  the  device  conductance  and  not  the 
susceptance. 

Figures  2.9  and  2.10  show  the  TE  device  admittance,  EF  voltage  and 
capacitance  as  a  function  of  frequency  for  optimum  load  (with  tuner)  and 
for  two  bias  values.  The  device  capacitance  is  calculated  from  susceptance 
and  is  seen  to  decrease  with  bias  increase  and  shows  a  tendency  to  increase 
with  frequency.  The  EF  voltages  are  larger  than  for  the  matched  load  case 
(Fig.  2.7)  but  still  are  insufficient  for  Q-mode  over  a  portion  of  X-band 
in  both  cases. 

Figure  2.11  shows  the  TE  device  characteristics  near  9  GHz  for 
various  values  of  bias.  A  decrease  in  conductance  and  susceptance  clearly 
occurs  for  increasing  bias  values.  Furthermore  EF  voltage  is  sufficient  for 
Q-mode  at  each  value  of  bias. 


FIG. 


.8  LOAD  CONDUCTANCE  AND  SUSCEPTANCE  VS.  FREQUENCY  FOR  TE  DEVICE 


M61-24  FOR  THE  MATCHED  LOAD  CASE  AND  FOR  THE  BEST  RF  POWER 
CONDITION  USING  AN  RF  TUNER  IN  LOAD  LINE.  (VB  =  15  V) 


DEVICE  CONDUCTANCE  CD,  pF  PEAK  RF 

AND  SUSCEPTANCE,  mm  ho  VOLTAGE, 


8  9  10  II  12  13 

FREQUENCY,  GHz 


FIG.  2.9  ADMITTANCE,  CAPACITANCE  AND  RF  VOLTAGE  AS  A  FUNCTION  OF  RF 

FREQUENCY  FOR  TE  DEVICE  M61-24  IN  THE  WAVEGUIDE  CIRCUIT  WITH 
A  TUNER  FOR  BEST  OUTPUT  POWER.  (Vn  =  15  V} 

'  iJ 
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FIG. 


2.10  ADMITTANCE,  CAPACITANCE  AND  RF  VOLTAGE  AS  A  FUNCTION  OF  RF 

FREQUENCY  FOR  TE  DEVICE  M6l-24  in  THE  WAVEGUIDE  CIRCUIT  WITH 
A  TUNER  FOR  BEST  OUTPUT  POWER.  (V-  =  20  V) 

D 
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FIG.  2.11  RF  POWER  OUTPUT,  FREQUENCY  OF  OSCILLATION,  RF  VOLTAGE,  LOAD 
CONDUCTANCE  AND  SUSCEFTANCE  AS  A  FUNCTION  OF  BIAS  VOLTAGE  FOR 
TE  DEVICE  Mol-24  IN  THE  WAVEGUIDE  CIRCUIT  WITHOUT  A  TUNER  AND 
FOR  A  POST-TO-SHORT  DISTANCE  OF  2.354  cm. 
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Conclusions.  Analysis  of  the  post-type  WR-90  waveguide  circuit 
has  shown  that  for  post-to-short  plane  separations  close  to  one-half 
wavelength  only  one  circuit  resonance  occurs  in  X-band  for  capacitive  Gunn 
devices.  This  occurs  above  a  frequency  of  series  resonance  and  slightly 
below  a  frequency  of  parallel  resonance.  For  a  matched  load  condition,  the 
real  loading  of  the  TE  device  increases  greatly  for  an  increase  of  circuit 
resonant  frequency.  This-  is  responsible  for  the  observed  power  fall  off 
with  frequency  for  TE  device  m6i-24.  Use  of  a  slide-screw  tuner  was  shown 
to  result  in  lower  load  conductance  and  larger  EF  power  over  all  of  X-band. 

Careful  measurements  of  device  admittance  and  EF  voltage  verified 
the  preliminary  data  obtained  by  other  methods  and  reported  in  previous 
quarterly  progress  reports.  In  particular,  EF  voltage  is  found  to  be 
sufficient  for  Q-mode  operation  only  when  load  conductance  is  reduced  to 
values  nearly  optimum  for  best  EF  power  output.  It  is  possible  that  harmonic 
voltages  aid  the  quenching  process  for  the  higher  load  condition. 

The  operating  characteristics  at  15  V  bias  and  below  are  most 
consistent  with  Q-mode  behavior.  Operation  at  higher  voltages  probably 
involves  incomplete  domain  formation.  Characterization  of  the  active  diode 
as  a  parallel  capacitor  and  a  negative  resistor  is  appropriate  over  a  wide 
frequency  range.  For  fixed  voltage  bias,  the  capacitance  is  not  a  strong 
function  of  frequency  or  EF  voltage.  The  capacitance  reduces  as  bias  is 
increased  and  is  about  twice  the  low-field  value  at  20  V  which  is  ^.5  times 
threshold.  This  result  is  consistent  with  Q-mode  theory.  The  magnitude 
of  the  device  negative  conductance  decreases  with  increase  of  either  bias 
voltage  or  RF  voltage. 
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2.5  Program  for  the  Next  Quarter.  The  experimental  study  of  pulsed 
and  CW  TE  oscillators  will  be  continued.  The  circuit-device  interaction 
problem  wi  L  ae  studied  in  an  effort  to  understand  why  RF  efficiency  varies 
sc  much  as  a  function  of  frequency. 

2*  Gunn-Effect  Amplifiers 

Supervisor:  W.  R.  Curtice 
Staff:  A.  K.  Talwar 

j.l  Introduction.  In  Quarterly  Progress  Report  No.  3>  numerically 
calculated  results  of  a  large-signal  analysis  of  stabilized  transferred- 
electron  devices  were  given.  The  dependence  of  impedance  and  efficiency 
upon  RF  signal  amplitude  and  length  of  the  device,  etc.,  was  described. 
Effects  of  diffusion  were  also  discussed.  In  this  report,  plots  of  small- 
and  large-signal  impedance  against  frequency  are  shown  in  order  to  enable 
an  assessment  of  the  bandwidth  capabilities  of  these  devices.  Nearly  octave 
bandwidths  as  observed  experimentally  by  Perlman"  have  been,  calculated.  It 
is  also  shown  in  this  report  that  inclusion  of  field  dependence  of  diffusion 
is  essential  for  an  accurate  evaluation  of  the  dc  characteristics  of  these 
devices.  The  effects  of  lattice  temperature  on  the  various  characteristics 
of  these  devices  is  investigated.  It  is  shown  that  a  rise  in  lattice 
temperature  does  not  necessarily  reduce  the  magnitude  of  negative 
resistance.  Sufficient  increase  in  lattice  temperature  above  the  room 
temperature  may  be  necessary  for  significant  reduction  of  the  magnitude  of 
negative  RF  resistance  in  order  to  aid  stabilization  of  the  devices.  The 

1.  Perlman,  B.  S.,  "Microwave  Amplification  Using  Transferred-Electron 
Devices  in  Prototype  Filter  Equalization  Networks,"  RCA  Rev.,  vol.  32, 
No.  1,  pp.  3-23,  March  1971- 
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effects  of  gradients  in  lattice  temperature  and  donor  density  along  the 
length  of  a  stabilized  GaAs  crystal  are  discussed  analytically.  It  is 
suggested  that  such  gradients  can  have  significant  influence  on  the 
performance  of  these  devices.  In  particular,  RF  to  dc  conversion 
efficiency  can  be  improved  by  properly  controlling  such  gradients. 

3.2  Bandwidth  Considerations,  figure  3*1  shows  the  resistance  and 
reactance  of  a  l'+-(jm  long  crystal  with  nQ  =  5  x  10l4/cm3  plotted  as  a 
function  of  frequency  for  different  values  of  normalized  RF  current.  As 
before,  the  currents  are  normalized  with  respect  to  as 


th 


n  e  v 
o 


th  ’ 


(3.1) 


where  v,,  is  the  carrier  velocity  at  the  threshold  electric  field,  J  is 
th  c 

taken  to  be  10  percent  above  and  the  area  was  taken  to  be  9*8  x  10  4  cm2. 

The  impedance  results  show  two  negative  resistance  frequency  regions 

which  may  be  used  for  reflection  amplification.  The  low-frequency  region 

occurs  near  the  transit-time  frequency  since  6  GHz  times  l4  urn  is 

0.84  x  107  cm/s.  Note  here  that  the  reactance  is  about  twice  the  value  of 

(negative)  resistance.  This  relationship  influences  the  gain-bandwidth 

product  when  the  device  is  used  for  amplification. 

J  =  0.001  is  a  small-signal  case,  since  any  value  less  than  this 
n 

leads  to  the  same  results.  It  is  seen  that  the  maximum  negative  resistance 

at  first  increases,  as  J  is  increased,  and  then  decreases.  Figure  3-2 

n 

shows  a  plot  similar  to  that  of  Fig.  3*1  for  a  crystal  with  the  same  donor 
density,  but  with  a  length  of  10  ,im  and  an  area  of  7  x  10  4  cm2.  The 
low-field  resistance  is  kept  the  same  (nearly  2.5  H)  in  the  two  cases. 

The  dc  current  density  is  also  the  same.  The  shorter  device  (Fig.  3*2) 
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mu 


has  a  slightly  smaller  magnitude  of  negative  resistance,  As  expected,  the 
negative  resistance  frequency  region  is  shifted  to  higher  frequencies.  The 
frequency  range  of  this  region  is  larger  than  that  for  the  longer  device. 
The  last  result  seems  to  agree  with  the  calculations  of  Magarshak  and 
Mircea2  who  have  found  very  large  bandwidths  for  short  crystals.  The 
shapes  of  the  resistance  and  reactance  plots  of  Figs.  3*1  and  3*2  bear  a 
marked  resemblance  to  the  experimental  curves  of  Perlman.1 

3.3  Field  Dependence  of  Diffusion  Coefficient  and  Its  Importance  in 
Determining  the  dc  Characteristics  Accurately.  The  calculated  value  of 
average  carrier  density  n,  and  the  values  of  diffusion  coefficient 
corresponding  to  the  average  electric  field  at  each  point  along  the  length 
of  the  crystal  of  Fig.  3-1  are  plotted  in  Fig.  3*3  for  a  small  RF  signal 
-evel.  Note  that  the  dc  current  due  to  a  gradient  in  carrier  density  and 
that  due  to  a  gradient  in  diffusion  coefficient  flow  in  opposite  directions 
Thus  the  assumption  of  a  field-independent,  constant  diffusion  coefficient 
would  lead  to  erroneous  results.  Clearly,  a  more  accurate  dc  solution  is 
obtained  if  the  diffusion  is  neglected  completely  than  if  it  is  taken  to  be 
a  constant.  Of  course,  taking  into  account  the  field  dependence  of  the 
diffusion  coefficient  gives  the  most  accurate  results. 

hi  Effects  of  Lattice  Temperature.  In  this  section  it  is  assumed 
that  the  lattice  temperature  is  the  same  everywhere  inside  the  crystal.  An 
increase  in  lattice  temperature  reduces  both  the  low-field  mobility  and 
the  negative  differential  mobility  and  changes  the  static  field  profile  in 

2.  Magarshak,  J.  and  Mircea,  A.,  "Wideband  CW  Amplification  in  X-Band 

with  Gunn  Diode',"  Int.  Solid  State  Circuits  Conf .  Digest,  Philadelphia 
Pa.,  p.  134,  February  1970* 
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the  crystal.  In  order  to  understand  the  effects  of  these  changes,  small-signal 
impedance  calculations3  were  made  assuming  a  piece-wise  linear  v-E  character¬ 
istic  and  D  =  0.  Figure  3  A  shows  the  different  v-E  relationships  used  and 

these  have  different  peak- to -valley  ratios.  Computations  were  made  for  a  crystal 

of  10  jim,  nQ  =  5  x  1014/cm3  and  an.-i  -  7  x  10  4  cm2  over  a  band  of  frequencies 

and  the  maximum  negative  i  distance  obtained  (around  10  GHz)  is  plotted 

against  peak-to -valley  ratios  ir.  Fig.  3*5  for  a  value  of  Jq  equal  to  1.1. 

Maximum  value  of  negative  resistance  in  the  frequency  band  is  chosen  so  that 

transit  angle  effects  are  unimportant  for  the  following  discussion.  Although 

a  direct  decrease  in  the  magnitude  of  negative  resistance  with  a  decrease  in 

the  j  k- to -valley  ratio  was  expected,  it  did  not  always  occur.  Instead, 

the  maximum  negative  resistance  was  not  greatly  affected  by  changes  in  the 

peak-to -valley  ratio  unless  the  ratio  was  small.  The  reason  for  this 

behavior  can  be  explained  with  the  aid  of  Fig.  3*5  which  shows  the  values 

of  L,  the  length  of  the  active  region,  and  L  ,  the  length  of  the  anode 

region.  A  large  peak-to -valley  ratio  results  in  a  large  blowing  down  of 

charge  carriers  with  a  consequent  large  static  field  distortion.  This  leads 

to  a  smaller  length  L  of  the  active  region  and  a  longer  passive  region  (L  ). 

n.  a. 

A  decrease  in  the  peak-to -valley  ratio  increases  the  length  of  the  active 

region  (L  )  and  decreases  the  length,  of  the  passive  region  (L  ).  This 
n  a 

compensates  for  the  reduced  negative  differential  mobility  until  lower  values 
of  peak-to -valley  ratio  occur.  Thus  a  significant  increase  in  lattice 
temperature  above  room  temperature  is  necessary  to  reduce  the  magnitude  of 
the  negative  resistance  by  means  of  a  change  in  the  v-E  curve. 

3.  Haddad,  G.  I.  et  al.,  "Microwave  Solid-State  Device  and  Circuit  Studies," 
Tech.  Report  RADC-TR-71-109,  Electron  Physics  Laboratory,  The  University 
of  Michigan,  Ann  Arbor,  June  1971* 
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FIG.  5.4  VELOCITY-FIELD  RELATIONSHIPS  USED  FOR  SMALL-SIGNAL  DIFFUSION- 
FREE  STUDY. 
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PEAK-TO-VALLEY  RATIO 


FIG.  3.5  DIFFUSION -FREE,  SMALL-SIGNAL  CALCULATIONS  OF  MAXIMUM  NEGATIVE 

RESISTANCE  AND  LENGTH  OF  ACTIVE  (L  )  AND  PASSIVE  (L&)  REGIONS  IN 

A  DEVICE  10  um  LONG  WITH  nQ  =  5  x  10l4/cm3  AND  AREA  =  7  x  10"4  cm2. 

VELOCITY-FIELD  RELATIONS  OF  FIG.  3.^  WERE  USED. 
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It  was  mentioned  previously  that  static  field  distortion  is  high 

when  the  peak-to -valley  ratio  is  large  or  at  lower  temperatures.  To 

illustrate  this  effect,  plots  of  static  electric  field  vs.  distance  from  the 

cathode  are  shown  in  Fig.  3*6  for  n  =  5  x  10l4/cm3  and  J  =1.1.  These 

o  o 

were  calculated  by  applying  the  large -signal  analysis  including  diffusion 
effects,  for  the  analytical  forms  of  v-E  characteristics,  described  in  the 
last  quarterly  progress  report  for  300°K,  375°K  and  500°K.  These  calcu¬ 
lations  are  for  a  small-signal  case. 

The  large-signal  resistance  and  reactance  were  then  calculated  for 
the  same  material  parameters  and  a  bias  of  approximately  l6  V:  l  =  10.1  qm, 
area  =  7  x  10  4  cm2  and  J  ,  the  normalized  RF  current,  equals  0.15.  These 

li 

results  are  shown  in  Fig.  3*7*  The  bias  voltage  was  kept  constant  by 
adjusting  the  value  of  Jq  for  each  temperature.  The  increase  in  temperature 
is  seen  to  reduce  the  magnitude  of  maximum  negative  resistance  and  also 
the  frequency  at  which  it  occurs.  The  major  changes  occur  at  this  highest 
temperature,  500°K. 

It  was  shown  earlier3  that  an  increase  in  RF  signal  drive  increases 
the  static  field  distortion  in  the  crystal  leading  to  a  greater  voltage 
drop  across  the  device  for  a  constant  dc  current.  In  a  constant  voltage 
circuit,  then,  as  the  RF  signal  level  increases,  the  increase  in  static 
field  distortion  will  lead  to  a  current  drop  in  the  bias  current. 

In  addition  to  this  current  drop,  Perlman  et  al.4  have  observed  that 
in  CW  devices  additional  current  drop  occurs  at  high  bias  voltages  due  to 
heating  of  the  crystal.  The  present  theory  permits  a  calculation  of  this 

4.  Perlman,  B.  S.,  Upadhyayula,  C.  L.  and  Siekanowicz,  W.  W.,  "Microwave 
Properties  and  Application  of  Negative  Conductance  Transferred-Electron 
Devices,"  Proc.  IEEE,  voi.  59>  No.  8,  pp.  1229-1237,  August  1971* 
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DISTANCE  FROM  CATHODE,  /im 

FIG.  3.6  STATIC  ELECTRIC  FIELD  VS.  DISTANCE  FROM  THE  CATHODE  FOR  A  SMALL- 
SIGNAL  CASE  AT  DIFFERENT  LATTICE  TEMPERATURES  WITH  DIFFUSION 
EFFECTS  INCLUDED.  (Jq  =  1.1,  J  =  0.001,  FREQUENCY  =  10  GHz  AND 
nQ  =  5  x  1014/cm3) 
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current  drop  and  some  results  are  shown  in  Fig.  3.8  where  the  device  dc 
current  is  plotted  as  a  function  of  bias  voltage  for  small  ac  signals  and  a 
crystal  length  of  10  |im,  nQ  =  5  x  10l4/cm3  and  area  =  3*5  x  10  4  cm2.  The 
plot  is  in  two  parts.  The  section  at  low  bias  voltages  is  for  a  lattice 
temperature  of  375°K  and  the  one  at  higher  bias  voltages  is  for  a  lattice 
temperature  of  500°K.  If,  in  the  CW  operation,  the  bias  is  increased 
resulting  in  an  increase  in  the  lattice  temperature  from  375 °K  to  500°K, 
the  operation  of  this  device  will  move  from  some  point  on  the  375°K  curve 
to  a  point  on  the  500°K  curve  depending  on  the  thermal  resistance  of  the 
packaged  device.  If  the  device  is  initially  operated  at  point  A  in  the 
figure,  then  the  bias  is  increased  and  the  lattice  temperature  reaches  500°K. 
The  operating  points  for  four  different  values  of  the  thermal  resistance  are 
shown  in  the  figure.  All  show  a  current  drop  from  point  A.  The  higher  the 
thermal  resistance,  the  smaller  is  the  bias  voltage  required  to  reach  500°K 
and  the  larger  is  the  current  drop. 

hi  Variation  of  the  Lattice  Temperature  Along  the  Length  of  the 
GaAs  Device .  It  was  shown  in  the  last  quarterly  progress  report  that  an 
increase  in  static  field  distortion  causes  an  increase  in  dc  power 
dissipation  and  consequently  a  reduction  in  the  dc  to  RF  conversion 
efficiency  of  stabilized  GaAs  devices.  The  static  field  distortion  occurs 
in  GaAs  because  the  average  carrier  velocity  is  smaller  at  high  electric 
fields  and  since  the  electric  field  increases  monotonically  from  cathode 
to  anode,  the  charge  carriers  have  a  small  average  velocity  near  the  anode 
compared  to  the  point  where  the  electric  field  is  near  threshold.  This 
difference  in  velocity  causes  charge  accumulation  to  occur  near  the  anode 
in  order  to  satisfy  the  condition  of  continuity  of  current.  If  the  region 
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FIG.  3.8  Be  DEVICE  CURRENT  VS.  BIAS  VOLTAGE  AT  TWO  DIFFERENT  LATTICE 
TEMPERATURES.  (LENGTH  =  10  pm,  AREA  =  3.5  x  10'4  cm2  AND 
nQ  =  5  x  1014/cm3) 
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of  the  crystal  near  the  anode  were  cooled  to  a  lower  temperature  in 
c= -^orison  to  the  region  near  the  cathode,  the  velocity  of  the  charge 
carriers  near  the  anode  would  be  increased  and  this  would  lead  to  a  reduced 
space-charge  accumulation  and  consequently  a  smaller  static  field  distortion. 
The  same  dc  result  can  be  obtained  by  heating  the  whole  crystal,  thereby 
reducing  the  peak -to -valley  ratio.  The  situation  is  not  the  same  however 
when  EF  behavior  is  considered.  In  the  dc  case  the  field  varies  only  in 
space,  whereas  when  RF  is  considered  there  is  also  a  fluctuation  of  electric 
field  with  respect  to  time.  Thus  it  is  likely  that  different  effects  are 
obtained  when  there  ire  lattice  temperature  gradients  with  temperature 
decreasing  away  from  the  cathode  than  when  the  whole  crystal  lattice  is 
heated. 

A  study3  of  the  effects  of  lattice  temperature  can  be  made  by  assuming 

the  parameters  a  ,  a  and  E  of  the  velocity- field  characteristics  to  be 
\  2  n 

dependent  on  x,  the  distance  from  the  cathode.  Preliminary  calculations 
have  been  made  and  it  has  bean  found  that  considerable  improvement  in  the 
efficiency  can  be  achieved  if  the  lattice  temperature  decreases  away  from 
the  cathode.  The  results  will  be  reported  in  a  future  report. 

3.6  Variation  in  Donor  Density  Along  the  Length  of  a  GaAs  Crystal. 

It  has  been  stated  previously  that  the  static  field  distortion  results  from 
an  accumulation  of  charges  near  the  anode,  required  by  the  negative  differ¬ 
ential  mobility  in  GaAs.  It  was  maintained  in  Section  3*5  that  the 
accumulation  can  be  reduced  by  creating  a  gradient  in  lattice  temperature 
along  the  length  of  the  device.  Another  way  of  reducing  the  accumulation  of 
excess  carriers  would  be  to  increase  the  donor  aensity  away  from  the  cathode. 
It  is  thus  suggested  that  efficiency  should  improve  if  the  doping  profile 
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of  the  GaAs  crystal  is  properly  controlled.  In  the  analysis  presented  in  the 

previous  quarterly  progress  reports,  the  donor  density  nQ  was  assumed  constant. 

If  n  is  not  constant,  then  the  diffusion  current  density  must  he  written  as 
o 


-e 


dE  \  dD  3e 
<5x  )  ^E  ctoc 


dx2 


3a 
_ o 

ST  * 


(3-2) 


The  last  term  is  the  new  contribution.  The  solution  of  the  problem  with 
3nQ/3x  ^  0  can  be  carried  out  as  done  previously  in  Quarterly  Progress 
Report  No.  3* 

3 »7  Conclusions.  Consideration  of  many  numerical  cases  using  the 
large-signal  theory  of  stabilized  TE  devices  permits  the  following  conclusions 
to  be  drawn: 

1.  Stabilized  TE  devices  show  a  fairly  large  bandwidth  which  appears 
to  increase  as  the  length  of  the  crystal  is  reduced. 

2.  Inclusion  of  the  field  dependence  of  diffusion  is  important  for 
accurate  evaluation  of  dc  quantities. 

3*  High  operating  temperatures  are  necessary  for  significant  reduction 
of  the  magnitude  of  negative  RF  resistance  which  aids  stabilization  of  the 
TE  devices. 

4.  it  is  suggested  that  static  field  distortion  in  low  nQi  devices 
can  be  reduced  by  creating  a  lattice  temperatura  gradient  with  temperature 
decreasing  away  from  the  cathode  or  by  increasing  the  doping  density  near 
the  anode. 

3.8  Program  for  the  Next  Quarter.  Theoretical  calculations  will 
be  continued  to  study  the  effect  of  varying  donor  density  in  a  stabilized 


-33- 


GaAs  transferred-electron  device.  Experiments  will  be  conducted  to  verify 
the  salient  results  of  the  large-signal  analysis.  Experiments  in  the 
delayed-domain  mode  of  operation  will  also  be  conducted. 

4 .  Modulation  Properties  of  Gunn-Effect  Devices 
Supervisor:  W.  R.  Curtice 
Staff:  D.  Tang 

4.1  Introduction.  The  objective  of  the  present  study  is  to 
determine  the  effect  of  lattice  temperature  upon  the  frequency  modulation 
sensitivity  of  Gvnn-eifect  devices.  The  experimental  results  are  presented 
and  discussed  in  the  following  section. 

In  addition,  the  impedance  of  a  waveguide  mounting  structure  is 
investigated  for  the  purpose  of  designing  a  cavity  that  has  constant  real 
loading  over  X-band  frequency  and  a  large  tuning  rate. 

4.2  Experimental  Investigation  of  Modulation  Properties  of  a  Gunn- 
Effect  Device  at  Different  Temperatures.  Experiments  at  different 
temperatures  were  performed  using  the  Gunn-effect  device  HB-10B-105  mounted 
in  a  50-fi  coaxial-line  cavity.  The  cavity  slugs  are  set  so  chat  resonant 
frequency  of  the  cavity  is  about  8.6  GHz  and  the  slugs  are  not  removed 
throughout  the  experiments.  The  experimental  setup  is  drawn  schematically 
in  Fig.  4.1, 

The  I-V  characteristics  are  measured  ana  plotted  in  Figs.  4.2  through 
4.4  for  heat  sink  temperatures  of  8.5°c,  -38°C  and  -6l°C,  respectively.  The 
threshold  current  increases  and  low-field  resistance  decreases  as  temperature 
drops.  This  reflects  the  mobility  increases  when  temperature  drops  because 
the  electron  carrier  concentration  is  constant  through  the  temperature  range 
of  the  experiments.  The  large  amount  of  dc  current  decrease  at  large  bias 
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BIAS  CIRCUIT 


FIG.  4.1  EXPERIMENTAL  SETUP. 


BIAS  CURRE 


FIG.  4.2  I-V  CHARACTERISTICS  OF  DEVICE  HB-10B-105  FOR  HEAT  SINK  TEMPERATURE 
AT  8.5°C.  THE  LOW-FIELD  RESISTANCE  IS  8  fl. 


-56- 


L 


voltage  results  from  the  device  large  thermal  resistance.  The  device 
temperature  increases  rapidly  with  bias  increase  resulting  in  a  rapid 
decrease  in  current  with  an  increase  in  bias  voltage.  Device  parameters 
are  estimated  to  be  as  follows: 
device  length  -  10  urn, 

carrier  concentrations  1  to  3  x  10lS  cm  3, 

Hall  mobility  at  room  temperature  -  8.2  x  103  cm2/v-s  and 
thermal  resistance  -  26o°c/w. 

Data  for  frequency  modulation  sensitivity  (FMS)  is  taken  at  different 
bias  voltages  for  temperatures  at  8.5°C,  -2l°C,  -38°C  and  -6l°C.  The  data 
show  the  FMS  is  constant  over  the  modulation  frequency  range  from  200  kHz 
to  800  kHz.  It  is  also  constant  over  the  modulation  voltage  from  0  V  to 
1.5  V  (peak  value).  Typical  data  are  shown  in  Figs.  4.5a  and  4.5b.  The  FMS 
varies  rapidly  with  temperature  and  bias  voltage.  Its  value  varies  from 
1  MHz/v  to  17.6  MHz/v.  Figure  4.6  is  a  plot  of  FMS  vs.  bias  voltage  for 
different  heat  sink  temperatures.  Figure  4.7  is  a  plot  of  FMS  vs.  device 
average  temperature  by  assuming  thermal  resistance  of  the  device  to  be 
260°c/w.  Generally,  FMS  decreases  w  en  either  temperature  or  bias  voltage 
is  increased. 

Amplitude  modulation  sensitivity  (AMS)  is  defined  as  the  percentage 
power  variation  per  unit  modulation  voltage  (peak  value).  All  the  data 
show  that  AMS  decreases  as  modulation  frequency  increases  as  shown  in 
Figs.  4.8  and  4.9.  The  amplitude  envelope  of  AM  waveform  is  observed 
from  the  CRO.  It  is  found  that  at  low  modulation  frequency  the  waveform 
is  more  distorted  than  at  high  modulation  frequency.  This  occurs  because 
the  device  is  nonlinear  with  respect  to  the  modulation  voltage.  Higher 
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PIG.  4.5  FREQUENCY  MODULATION  SENSITIVITY  FOR  HEAT  SINK  TEMPERATURE  AT  -38°C 
(a)  AS  A  FUNCTION  OF  PEAK  MODULATION  VOLTAGE  AND  (b)  AS  A  FUNCTION 
OF  MODULATION  FREQUENCY. 
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FREQUENCY  MODULATION  SENSITIVITY,  MHz/V 
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FIG.  4.7  FREQUENCY  MODULATION  SENSITIVITY  VS.  DEVICE  AVERAGE  TEMPERATURE. 
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FIG.  fc.8  AMPLITUDE  MODULATION  SENSITIVITY  FOR  HEAT  SINK  TEMPERATURE  AT  8.S°C 
AND  RF  FREQUENCY  AT  8.56  GHz.  (a)  AS  A  FUNCTION  OF  MODULATION 
FREQUENCY  WHEN  THE  MODULATION  VOLTAGE  IS  AT  0.2  V  AND  fb')  AS  A 
FUNCTION  OF  MODULATION  VOLTAGE. 


(b) 

FIG.  4.9  AMPLITUDE  MODULATION  SENSITIVITY  FOR  HEAT  SINK  TEMPERATURE  AT  -6l°C 
AND  RF  FREQUENCY  =8.68  GHz.  (a)  FOR  PEAK  MODULATION  VOLTAGE  AT 
0.6  V  AS  A  FUNCTION  OF  MODULATION  FREQUENCY  AND  (b)  AS  A  FUNCTION 
OF  PEAK  MODULATION  VOLTAGE  AT  DIFFERENT  BIAS  VOLTAGES. 


harmonics  of  the  modulation  frequency  are  produced  by  the  device.  The 
distortion  is  more  apparent  at  low  temperature. 

4.3  Device  Operation  in  Waveguide  and  Coaxial  Circuits.  In  this 
section  the  modulation  sensitivities  are  discussed  qualitatively.  By 
plotting  the  cavity  circuit  admittance  and  the  negative  of  device  admittance 
on  the  complex  admittance  plane  it  can  be  seen  how  the  circuit  and  the 
device  affect  the  modulation  sensitivities.  The  coaxial  cavity  and  the 
waveguide  cavity  circuit  have  some  similarities. 

Figure  4.10a  is  a  plot  of  the  circuit  line  of  the  coaxial  cavity  used 
in  the  experiment  and  Fig.  4.11a  is  a  typical  X-band  waveguide  cavity 
circuit  line.  Both  lines  are  computed  including  the  package  reactances. 

The  scales  on  the  circuit  lines  mark  the  corresponding  frequency.  The 
scale  is  wider  at  the  right-hand  part  of  the  circuit  line.  The  negative  of 
the  device  line  is  also  drawn  in  Figs.  4.10a  and  4.11a.  The  stable 
oscillation  point  is  the  intersection  of  the  device  line  and  the  circuit  line 
at  which  the  circuit  line  satisfies  the  condition  £©/df  >  0.  In  Figs.  4.10a 
and  4.11a  the  points  labeled  2  are  the  stable  operating  points.  Figures 
4.10b  and  4.11b  are  the  enlarged  plot  of  the  region  near  the  operating  point 
in  Figs.  4.10a  and  4.11a,  respectively. 

From  the  experimental  data  for  the  coaxial  circuit  it  is  found  that 
increasing  the  bias  voltage  usually  causes  the  frequency  to  increase.  This 
reflects  that  the  intersection  point  shifts  upward  along  the  circuit  line. 

If  the  modulation  frequency  is  not  so  high  that  the  device  admittance  behaves 
as  it  does  at  steady  state,  it  is  found  that  the  operating  point  traces  up 
and  down  along  a  section  of  the  circuit  line.  This  is  shown  in  Fig.  4.11b. 
The  FMS  is  proportional  to  the  frequency  range  of  the  circuit  line  covered 
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FIG.  4.10  (a)  CALCULATED  COAXIAL  CAVITY  CIRCUIT  AND  DEVICE  LINES  FOR  THE 

EXPERIMENT  AND  (b)  ENLARGED  FIGURE  NEAR  POINT  2  IN  FIG.  4.10a. 
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4.11  (a)  TYPICAL  WAVEGUIDE  CAVITY  ADMITTANCE  AND  NEGATIVE  DEVICE 

ADMITTANCE  FCR  A  CONSTANT  BIAS  VOLTAGE  AND  (b)  ENLARGED  FIGURE 
NEAR  POINT  2  IN  FIG.  4.11a  AND  DEVICE  LINES  FOR  DIFFERENT  BIAS 


VOLTAGES. 


by  that  section.  The  AMS  is  proportional  to  tan  0,  where  9  is  the  intersect 
angle  between  the  constant  RF  voltage  contour  on  the  device  lines  and  the 
circuit  line. 

Employing  the  preceding  admittance  plot,  it  can  also  be  seen  how  the 
device  temperature  affects  the  modulation  sensitivity.  Raising  the  device 
temperature  always  decreases  the  oscillation  frequency.  This  means  the 
operating  point  shifts  downward  along  the  circuit  line.  If  the  scale  is 
wider  there,  smaller  FMS  is  expected. 

Increasing  bias  voltage  is  usually  accompanied  by  raising  device 
temperature  because  the  dc  power  dissipation  increases.  The  resultant  effect 
of  increasing  bias  may  shift  the  operating  point  either  upward  or  downward 
depending  on  the  thermal  resistance  of  the  device. 

Device  effects  can  be  discussed  by  comparing  the  experimental  data 
of  FMS  with  the  circuit  line  in  Fig.  4.10b.  Table  4.1  contains  the  data 
taken  for  heat  sink  temperature  at  -21 °C.  The  RF  frequency  changes  less 
than  70  MHz.  In  this  frequency  range,  the  frequency  scale  on  the  circuit 
line  is  almost  uniform  as  shown  in  Fig.  4.10b.  But  the  FMS  changes  nearly 
four  times.  It  is  apparent  that  the  difference  comes  from  the  device  itself. 


Table  4.1 

Gunn  Device  Behavior  for  Several  Bias  Voltages 


Bias  Voltage  (V) 

Fin 

RF  Frequency  (GHz) 

RF  Power 

7.91 

17-5 

8.614 

10.1 

8.4l 

13.0 

8.620 

12.5 

9-31 

4.4 

8.554 

13.0 
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4.4  Waveguide  Mounting  of  Devices.  In  a  waveguide  cavity  the  mounting 
post  works  as  an  antenna  in  the  waveguide.  KF  power  generated  from  the 
device  is  radiated  through  the  post  into  the  waveguide.  The  device  driving 
point  impedance  is  equal  to  the  radiation  impedance  of  the  post.  It  is 
found  that  post  reactance  is  directly  proportional  to  the  waveguide  height 
and  decreases  with  increasing  post  diameter.1  For  the  post-gap  structure 
which  is  used  for  mounting  the  device,  only  limited  references  are  available. 
Eisenhart's  theory2  is  applicable  to  the  case  where  the  gap  is  smaller  than 
one  fourth  of  the  guide  height  (his  post-gap  reactance  deviates  from  that 
given  by  Marcuvitz  by  22  H  for  the  gap  being  one  fourth  of  the  guide  height). 
When  the  device  and  its  package  are  put  into  the  gap,  the  strength  and  the 
distribution  of  the  electric  field  in  the  gap  are  different  from  that  of  the 
air  gap,  and  therefore  the  driving  point  impedance  is  different.  The  effect 
of  the  package  and  the  device  in  the  gap  changes  the  gap  capacitance.  It  is 
expected  that  if  the  package  height,  which  is  the  gap  height,  is  small 
compared  with  the  guide  height,  then  the  reactance  of  the  full  height  post 
is  an  approximation  for  the  mounting  post  and  the  gap  capacitance  must  be 
modified  bo  be  a  larger  value.  Its  exact  numerical  value  depends  on  the 
package  dimension. 

In  order  to  determine  the  circuit  parameters  which  limit  the  r/q 
ratio  (or  tuning  rate)  of  half-wavelength-type  waveguide  cavity,  calculations 


1.  Marcuvitz,  N.,  Waveguide  Handbook,  McGraw-Hill  Book  Co.,  Inc.,  New 
York,  pp.  271-272,  1951. 

2.  Eisenhart,  R.  L.,  "Theoretical  and  Experimental  Analysis  of  a  Waveguide 
Mounting  Structure,"  Ph.D.  Thesis,  The  University  of  Michigan,  1971* 
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were  made  with  circuit  parameters  changed  arbitrarily  and  the  R/Q  ratios 
were  computed  for  each  set  of  the  parameters.  It  was  found  that  neither 
the  post  inductance  nor  the  gap  capacitance  limits  the  r/q.  The  real  loading 
and  the  package  reactance  do  affect  r/q.  The  heavily  loaded  cavity  (lower 
values  of  characteristic  irapeaance  ZQ)  has  a  narrower  frequency  scale  and 
therefore  a  larger  tuning  rate.  For  a  reduced  height  X-band  waveguide 
cavity  with  height- to -width  ratio  of  ^.5,  R/Q  is  3.8  fi,  while  for  a  WR-90 
waveguide  cavity,  R/Q  is  about  unity.  The  tuning  rate  was  computed  for 
direct  coupling  of  the  package  into  a  fully  reduced  height  waveguide  cavity 
which  does  not  have  a  post  in  it,  and  the  waveguide  characteristic  impedance 
was  kept  the  same  by  reducing  the  width  of  the  waveguide.  The  calculated 
r/q  was  still  about  the  same  as  that  obtained  from  the  full  height  waveguide 
post  coupling  cavity.  Therefore  the  conclusion  is  reached  that  the  tuning 
rate  is  mainly  determined  by  the  package  parameters  and  the  waveguide 
impedance. 

^.3  Conclusions.  Tests  of  a  CW  Gunn-effect  device  show  tnat  FMS  is 
constant  with  respect  to  the  modulation  frequency  and  the  modulation  voltage. 
It  increases  as  temperature  and/or  bias  voltage  decreases.  However,  AMS 
decreases  as  the  modulation  frequency  increases.  A  distorted  power  waveform 
is  observed  at  lower  modulation  frequency.  The  nonlinear  relationship 
between  AMS  and  the  modulation  voltage  is  more  apparent  at  low.  r  temperature 
as  the  harmonic  modulation  sidebands  are  produced. 

Effect  of  bias  voltage  and  temperature  on  the  FMS  and  AMS  can  be 
interpreted  as  the  shifting  of  the  device  operating  point  along  the  circuit 
line  where  the  distance  between  the  frequency  scale  and  the  circuit  line 
slope  are  different.  The  experimental  data  show  large  FMS  variations  in  a 
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small  RF  frequency  range  where  the  frequency  scale  of  the  circuit  line 
is  nearly  uniform.  This  FMS  variation  comes  from  the  device  itself. 

When  the  package  is  postmounted  in  a  waveguide  cavity,  the  post 
inductance  is  nearly  the  same  as  a  full  height  post  in  a  waveguide.  Neither 
the  post  inductance  nor  the  gap  (mounting)  capacitance  limits  the  r/q  ratio 
of  a  waveguide  cavity.  The  waveguide  characteristic  impedance  and  the 
package  parameters  limit  the  r/Q  ratio.  A  reduced  height  cavity  heavily 
loads  the  device  and  improves  the  r/q,  ratio. 

J+.6  Program  for  the  Next  Quarter.  The  investigation  of  modulation 
properties  of  Gunn-effect  devices  will  be  continueu.  Further  tests  of 
the  temperature  effects  and  bias  voltage  effects  will  be  made  for  the 
waveguide  mounting  structure:  mechanics],  tuning  saturation  will  be  studied. 

A  cavity  design  with  separately  controllable  external  Q  and  tuning  rate 
r/q  will  be  sought. 

_5 .  Avalanche -Diode  Amplifiers 
Supervisor:  G.  I.  Haddad 
Staff:  R.  W.  La ton 

5«1  Introduction.  The  objective  of  this  phase  of  the  program  is 
to  s  udy  the  stability,  large -signal  gam  bandwidth  and  saturation  properties 
of  ivalancbe-diode  amplifiers  and  the  dependence  of  these  properties  upon 
doping  profile,  bias  mg  and  circuit  configuration. 

During  this  r  'porting  period  a  computer-aided  synthesis  procedure  was 
developed  for  the  determination  of  the  equivalent  circuit  characterizing  any 
particular  diode  package  and  mount  configuration.  Results  were  obtained 
utilizing  this  procedure  for  both  microstrip  and  ridged  waveguide  mounts. 


Wafer  admittance  characteristics  were  refined  using  the  Read  model  program 
described1  previously,  but  with  more  realistic  avalanche  region  widths  as 
determined  by  Schroeder  and  Haddad2  for  both  Si  and  GaAs  diodes.  Terminal 
plane  immittance  plots  were  obtained  using  the  transformation  determined  for 
the  packaged  diode  terminating  a  50-n  microstrip  line,  and  this  information 
was  used  to  begin  the  design  of  a  microstrip  reflection  amplifier  circuit. 

5.2  Determination  of  Immittance  Transformations  Characterizing  a 
Diode  Package  and  Mount.  In  the  previous  quarterly  progress  report  a  computer 
piugram  was  described  wnich  perturbs  assumed  equivalent  circuit  element 
values  in  such  a  manner  that  the  mean-square  error  between  measured  and 
calculated  terminal  impedance  values  is  minimized.  During  this  reporting 
period  a  systematic  method  of  determining  the  configuration  and  starting 
values  for  the  assumed  elements  has  been  developed,  based  upon  classical 
methods  for  the  synthesis  of  a  lossless  driving  point  impedance  from  the 
location  of  its  poles  and  zeros.3 * 

Begin  by  measuring  the  impedance  as  a  function  of  frequency  at  the 

desired  reference  plane  with  the  diode  biased  Ji.st  below  breakdown. 

Figure  5*1  is  a  plot  of  such  data  for  the  case  of  a  diode  in  a  Type  025 

package  terminating  a  0.650-inch  length  of  50-ft  microst'  ip  line  on  a  0.25- 

inch  thick  alumina  (e  -  10)  substrate.  From  this  figure  a  lossy  pole  of 

r 

1.  Schroeder,  W.  E.  and  Haddad,  G.  I.,  "Effect  of  Harmonic  and  Subharmonic 
Signals  on  Avalanche -Diode  Oscillator  Performance, "  IEEE  Trans .  on 
Microwave  Theory  and  Techniques  (Correspondence),  vol.  MTT-lB,  No.  6, 

PP.  327-331,  June  1970. 

2.  Schroeder,  W.  E.  and  Haddad,  G.  I.,  "Avalanche  Region  Width  in  Various 
Structures  of  IMPATT  Diodes,"  Proc.  IEEE  (Correspondence),  vol.  59, 

No.  8,  pp.  124-5-1248,  August  1QV1. 

3.  Guillemin,  E.  A.,  Synthesis  of  Passive  Networks,  John  Wiley  and  Sons, 

Inc.,  New  York,  Chap.  3,  1957* 
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FIG.  5.1  MEASURED  IMPEDANCE  FOR  A  DIODE  IN  AN  023  PACKAGE  TERMINATING  A 
50-fi  MICROSTRIP  LINE.  INSET  SHOWS  POSTULATED  POLE-ZERO 
CHARACTERIZATION. 
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impedance  is  clearly  evident  at  approximately  11.5  GHz.  A  passive  lossless 
circuit  must  have  poles  and  zeros  which  alternate  along  the  jo)  axis.  Thus 
the  simplest  realization  of  the  behavior  siiown  in  Fig.  5«1  is  obtained  from 
an  impedance  function  containing  a  zero  at  the  origin,  a  pole  at  11.5  GHz 
and  a  zero  at  infir ity.  In  fact,  however,  additional  data  at  frequencies 
above  and  below  those  plotted  in  Fig.  5*1  indicate  a  change  in  sign  of  the 
reactance  at  frequencies  immediately  below  7*5  GHz,  indicating  a  reactance 
zero  at  that  frequency  as  well  as  a  reactance  which  again  increases  with 
frequency  from  12  to  12.4  GHz.  The  pole-zero  configuration  is  postulated  as 
shown  at  the  top  of  Fig.  5*1.  The  impedance  function  may  then  be  constructed 
as 


Z(s) 


H 


(s2 


+  O)2  )(S2  +  CD2  ) 
Ol  _ 03__ 

s(r2  +  rf) 


(5-1) 


where  the  imredance  level  factor  H  is  chosen  to  match  reactances  at  some 
desired  s^  -  jox  near  the  center  of  the  frequency  range  of  interest. 

Expansion  of  Eq.  5*1  into  a  ladder  of  two  series  of  inductances  and  two 
shunt  capacitances  is  straightforward.  The  final  shunt  capacitance  should 
be  as  large  or  larger  than  the  independently  measured  Junction  capacitance. 

If  it  is  not,  the  impedance  level  of  the  final  loop  of  the  ladder  network 
may  be  adjusted,  for  example,  by  multiplying  the  appropriate  row  and  column 
of  the  impedance  matrix  by  a  constant.  Finally,  a  small  resistance  is  assumed 
in  series  with  the  junction  capacitance  to  represent  any  undepleted  substrate 
and  ohmic  contact  losses.  One  ohm  is  a  reasonable  value.  The  end  result, 
which  forms  the  starting  point  for  the  optimization  algorithm  by  the 
computer,  is  then  shown  in  Fig.  5-2.  Note  that  the  parallel  combination  of 
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FIG.  5-2  ASSUMED  EQUIVALENT  CIRCUIT. 
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C  and  C .  is  equal  to  the  capacitance  value  obtained  for  the  final  element 

1  J 

in  the  ladder  reduction  of  Fig.  5-1.  Starting  and  final  element  values  are 
shown  in  Table  5«1  and  the  calculated  real  and  reactive  parts  are  shown  by 
the  dashed  lines  in  Fig.  5*1. 

Table  5*1 

Equivalent  Circuit  Element  Values 
Element  Assumed  Starting  Value  Final  Value 


Rs 

1.5  ft 

1.5  ft 

c 

1 

0.188  pF 

0,2205  PF 

L 

2 

0.602  nH 

0.5^3  nH 

c 

2 

0.692  pF 

0.^924  pF 

L 

3 

0.285  nH 

0.5016  nH 

C.  (measured) 

0.55  PF 

0.55  PF 

Several  sources  of  inaccuracy  must  be  commented  upon.  First,  in  the 
microstrip  media  with  er  -  10,  the  wavelength  is  decreased  to  approximately 
4o  percent  of  its  free-space  value.  Therefore  at  10  GHz  ^microstrip  is 
approximately  0.480  inch.  The  outer  diameter  of  the  diode  package  ceramic 
is  0.085  inch,  which  is  about  two  thirds  of  a  quarter-wavelength.  Therefore, 
attempts  to  model  the  transformation  with  frequency-independent  lumped 
elements  can  only  be  expected  to  be  successful  over  fairly  narrow  ranges  of 
frequency.  Second,  the  impedance  was  measured  using  slotted-line  techniques 
with  an  open  circuit  at  the  physical  location  of  the  edge  of  the  diode 
package  ceramic  for  a  reference.  The  VSWR's  obtained  using  the  open  circuit 
were  rather  low,  particularly  above  11  GHz,  being  of  the  order  of  7  to  10. 
This  means  that  the  microstrip  line  itself  was  quite  lossy  in  that  frequency 
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range.  The  magnitude  of  the  real  part,  as  well  as  the  reactive  part,  near 
the  pole  at  11.5  GHz  is  very  strongly  influenced  by  this  loss,  and  both 
parts  can  be  duplicated  very  closely  from  11  to  12  GHz  by  allowing  the 
series  loss  to  go  as  high  as  10  SI.  The  usefulness  of  a  circuit  with  such 
a  high  loss  is  doubtful,  therefore  the  circuit  element  values  in  Table  5-1 
are  regarded  as  a  good  approximation  to  the  measured  impedances  from  8  to  11 
GHz,  even  though  they  produce  a  pole  (of  much  higher  impedance)  at  a 
frequency  of  13 • 5  GHz  instead  of  11.5  GHz. 

5« 3  Package  Transformation  in  a  Ridged  Waveguide .  The  techniques 
described  in  Section  5*2  have  been  applied  to  the  case  of  a  diode  in  an 
023  package  mounted  in  ridged  waveguide.  The  measured  and  calculated 
impedances  are  shown  in  Fig.  5*3  with  the  guide  dimensions  shown  in  the 
inset.  The  equivalent  circuit  is  shown  in  Fig.  5-^«  Since  the  diode  is 
mounted  in  a  transmission  configuration  with  the  far  end  of  the  guide 
terminated  in  a  matched  load,  the  50 -SI  parallel  resistance  across  the  package 
terminals  presents  Zq  the  characteristic  impedance  of  the  guide  at  approxi¬ 
mately  9  GHz.  Since  all  the  loss  could  be  placed  in  a  single  Rg  or  R^,  an 
ambiguity  is  introduced  into  the  optimization  scheme  by  the  choice  shown, 
with  an  attendant  loss  of  effectiveness  in  minimizing  real  part  error. 
Nevertheless,  the  result  is  considered  good  and  both  losses  are  retained  in 
the  interest  of  duplicating  the  physical  situation.  Finally,  all  the  imped¬ 
ances  are  taken  with  a  plane  transverse  to  the  axial  dimension  of  the  guide 
through  the  center  line  of  the  diode  as  a  reference. 

5.^  Admittance  Characteristics  of  Si  and  GaAs  Diodes .  Recent  work 
by  Schroeder  and  Haddad2  indicates  significant  differences  in  avalanche 
region  width  in  complementary  Si  diodes  because  of  the  differing  ionization 
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FIG.  50  MEASURED  AND  CALCULATED  IMPEDANCES  FOR  A  DIODE  IN  AN  023  PACKAGE 
MOUNTED  IN  RIDGED  WAVEGUIDE. 
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FIG-  5 


Element 


Rp  C=  z0(f)J 


Value 
1  ft 
50  ft 
0.586  nH 
0.84  nH 
0,24  pF 
0.94  pF 
0.57  PF 
0.55  pF 


4  EQUIVALENT  CIRCUIT  FOR  A  TYPE  02J  PACKAGE  IN  RIDGED  WAVEGUIDE 
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rates  for  holes  and  electrons.  Their  results  indicate  an  avalanche  region 
width  of  approximately  percent  in  the  commonly  utilized  p+nn+  Si  diodes. 

Since  our  earlier  theoretical  work  was  based  on  admittance  characteristics 
obtained  assuming  a  20  percent  avalaflche  width,  refinement  of  those  curves 
was  desirable.  This  was  accomplished  utilizing  the  Read  model  computer 
program  for  an  avalanche  width  of  2.7  ^m  and  a  drift  length  of  4  pm.  For 
comparison,  data  were  also  obtained  for  a  GaAs  diode  with  an  avalanche  width 
of  1  |im  and  a  drift  length  of  3  Mm*  The  bias  current  density  was  1000  A/cm2 
in  both  cases,  while  the  Si  diode  area  wat  made  larger  so  that  equal 
reactances  were  obtained  at  the  peak  value  of  small-signal  negative  conductance . 
The  results  are  shown  in  Fig.  5*5- 

Since  the  GaAs  diode  clearly  provides  larger  values  of  negative 
conductance,  its  characteristics  were  transformed  through  the  equivalent 
circuit  of  Fig.  5-2  to  provide  a  basis  for  design  of  a  microstrip  reflection 
amplifier.  The  transformed  characteristics,  on  an  impedance  basis,  are 
shown  in  Fig.  5-6. 

3.3  Program  for  the  Next  Quarter.  The  design  and  fabrication  of  a 
microstrip  reflection  amplifier  will  be  completed  and  experimental  measure¬ 
ments  undertaken  during  the  next  period.  In  addition,  theoretical  methods  of 
loading  and  combining  two  transmission  stages  of  amplification  in  ridged 
waveguide  will  be  investigated. 

6.  High-Efficicncy  Avalanche  Diodes 

Supervisors:  R.  J.  Lomax  and  G.  I.  Haddad 
Staff:  C.  M.  Lee 

6.1  Introduction.  The  objectives  of  this  phase  of  the  program  are 
to  investigate  numerically  and  experimentally  avalanche  diodes  with  a  view 
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FIG.  5-6  TERMINAL  IMPEDANCE  FOR  A  PACKAGED  (TYPE  025)  GaAs  DIODE  IN  A 
50-n  MICROSTRIP  LIRE  AS  A  FUNCTION  OF  RF  VOLTAGE  ON  THE  WAFER 
AND  FREQUENCY.  (J  =  1010  A/cm2) 


to  devising  techniques  to  increase  the  power  output  and  efficiency  of 
avalanche -diode  oscillators  and  amplifiers.  The  arm  of  the  numerical 
investigation  is  to  develop  a  mathematical  model  which  handles  distributed 
avalanche,  nonsaturated  drift  velocities  and  unequal  ionization  rates  for 
electrons  and  holes,  and  then  to  use  this  mathematical  model  to  investigat  •' 
the  operation  of  avalanche  diodes  of  various  constructions  under  various 
o  erating  conditions.  The  aim  of  the  experimental  investigation  is  to  test 
the  validity  of  the  theoretical  results. 

6.2  IMPATT  Mode  Computer  Program.  This  computer  program  has  been 
successfully  developed,  evaluated  and  used  extensively  in  the  investigation 
of  the  performance  of  complementary  diodes.  Its  operation  has  been  explained 
in  detail  in  Quarterly  Progress  Report  No.  3*  It  can  handle  distributed 
avalanche  and  unequal  ionization  rates  for  electro::5  and  holes,  however  it 
requires  that  the  drift  velocities  of  carriers  remain  saturated.  During  this 
quarter,  the  effect  of  mesh  size  on  the  calculated  admittance  has  been 
investigated.  Figure  6.1  is  a  sample  of  the  many  curves  obtained  in  this 
investigation  which  shows  the  convergence  of  the  solution  as  the  mesh  size 
becomes  smaller  and  smaller.  These  curves  will  help  to  determine  the  mesh 
size  required  for  a  given  accuracy.  Also  during  this  quarter  the  method  of 
calculating  the  admittance  has  been  changed  to  accelerate  the  convergence. 
Previously,  the  total  current  was  computed  from  the  following  equation 

JT  =  e  +  q  •  V  •  (P  +  N)  ,  (6.1) 

then  the  admittance  was  obtained  by  Fourier  analyzing  the  JT«  Since  Eq.  6.1 
involves  a  numerical  differentiation  (namely,  SE/6t)  which  always 
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amplifies  enors  in  E,  the  calculated  tended  to  have  large  errors. 
Furthermore,  because  the  conductance  G  usually  has  about  one -tenth  the 
magnitude  of  the  admittance  Y,  the  relative  error  in  G  was  about  ten  times 
larger  than  that  in  Y  or  in  All  these  effects  tended  to  make  the 
convergence  rate  slow.  In  the  new  method  the  total  current  is  separated 
ir.to  cwo  components:  the  displacement  current  and  the  induced  current 
L .  The  displacement  current  is  due  to  the  depletion  layer  capacitance, 
hence  its  associated  admittance  is  the  pure  susceptance  jB^  K.ich  can  be 
calculated  from 


jBd  =  jo£  =  jooe/w  .  (6.2) 

The  induced  current  J\  is  the  space  average  of  the  space-charge  current: 


Its  associated  conductance  G^  and  susceptance  B^  can  be  found  by  Fourier 
analyzing  <L(t).  Since  Eq.  6.2  djes  not  have  any  finite  dif^^.nce  error 
and  Eq.  6.5  involves  only  numerical  integration,  the  calculated  conductance 
(G  =  G^)  and  susceptance  (B  =  B^  +  Eh  )  are  much  less  affected  by  numerical 
errors. 

As  a  partial  check  for  this  program  and  for  general  use  (e.g.,  see 
Section  7)>  a  computer  program  has  been  written  to  calculate  the  small-signal 
admittance  of  an  avalanche  diode  of  arbitrary  doping  profile.  This  has  the 
advantage  of  being  inexpensive  to  run  and  to  some  extent  the  calculation  is 
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self-checking. 


6.3  TRAPATT  Mode  Program.  During  TRAPATT  mode  operation  the  electric 
field  in  the  avalanche  diode  may  become  vary  small  or  even  negative,  hence 
the  assumption  of  saturated  drift  velocities  made  in  the  IMPATT  mode  program 
is  no  longer  valid.  Previously,  the  carrier  curients  were  calculated  first 
from  the  following  equations: 


and 


JP(x,t)  =  JPsat  + 

/0X(G-  !)• dx 

(6>M 

JN(x,t)  =  JNsat  + 

fx  (a-^)  ■  ^  ■ 

(6.5) 

Next  the  carrier  densities  were  calculated  from 


and 


P(x, t)  =  JP(x, t)/VP 


N  (x,  t)  =  JN(x,t)/VN 


(6.6) 

(6.7) 


Difficulties  arise  when  the  velocities  become  zero  and  Eqs.  6.6  .and  6.7 
become  undefined.  One  way  to  avoid  these  difficulties  is  to  introduce  a 
diffusion  term.  Hence  to  consider  variable  velocities  it  is  helpful  to 
consider  diffusion  at  the  same  time,  and  the  continuity  equations  become: 


and 


d2P  c5P  SVp  5P  _ 

-D  —  +  V  +  r  +  TT  =  lx 
^x2  P  OX  OX  dt 


82N  6n  V  _ 

-D  -  -  V  r—  “  N  r  rr  =  [j 

^^2  n  ox  ox  dt 


(6.8) 


(6.9) 


An  iterative  method  is  used  to  solve  Eqs.  6.8,  6.9  and  Foisson's  equation.  A 
FORTRAN  program  that  will  implement  this  method  is  currently  being  wiitten. 
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6.4  Conclusions.  The  IMPATT  mode  computer  program  has  been  'well 


checked  out  and  used  extensively,  however,  its  application  is  limited  to 
the  IMPATT  mode  case  in  which  the  electric  field  does  not  usually  fall 
below  the  saturation  field.  Development  of  the  TRAPATT  mode  program  which 
will  handle  variable  velocities  is  being  carried  out. 

6.5  Program  for  the  Next  Quarter.  The  development  of  the  TRAPATT 
program  will  be  continued.  Once  a  correct  mathematical  model  is  achieved, 
systematic  theoretical  and  experimental  investigations  will  be  carried  out 
to  optimize  power  output  and  efficiency. 

7 •  Nonlinear  Operating  Characteristics  of  IMPATT  Diodes 
Supervisor:  G.  I.  Haddad 

Staff:  W.  E.  Schroeder 

7.1  Introduction.  In  the  last  quarterly  progress  report  a  large- 
signal  IMPATT  diode  analysis  which  utilizes  the  concept  of  separable 
ionization  and  drift  regions  was  described.  During  this  reporting  period 
extensive  numerical  results  were  obtained  from  the  computer  program.  Several 
cases  of  interest  are  described. 

"f* 

•  2_  Nonpunch-Through  Silicon  Diode.  Consider  the  case  of  an  n  pp 
diode  where  the  doping  level  and  width  of  the  p-region  are  such  that  the 
region  is  not  entirely  deplete^  when  the  diode  is  reverse  biased  into 
breakdown.  In  the  undepleted  epitaxial  region  the  electric  field  is  low 
and  the  current  flow  is  essentially  ohmic.  This  situation  is  frequently 
modeled  as  a  resistance  in  series  with  the  impedance  of  the  active,  i.e., 
depleted  region  of  the  diode.  This  approach  is  simple  and  relatively 
accurate  for  ^mall  RF  drive  levels;  however  for  '’urge  drive  levels  the 
depletion  width  is  modulated  and  she  field  in  the  undep] eted  region  may 
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exceed  that  required  for  constant  mobility,  thus  complicating  the  Picture. 

An  alternative  approach  is  to  consider  the  entire  p-region  to  be  the  active 
region  and  to  find  the  current  by  solving  the  carrier  transport  equations 
throughout  the  region. 

Using  the  IMPATT  diode  analysis  described  in  the  last  quarterly  progress 
+  + 

report  an  X-band  Si  n  pp  diode  which  breaks  down  prior  to  punch  through  was 
modeled.  The  p-region  was  taken  to  be  uniformly  doped  at  5  x  1015  cm  3  with 
a  total  width  of  9  pm,  composed  of  a  1-um  avalanche  region  and  an  8-|am  ''drift" 
region.  The  ionization  rates  of  electrons  and  holes  were  taken  equal  to  a 
single  effective  rate  which  is  appropriate  to  Si  at  room  temperatures: 

“eff  =  A  exP  G  |)  ’ 

where  A  =  1.315  x  107  cm  1  and  b  =  2.46  x  v/cm.  The  expression  for 
drift  velocity  as  a  function  of  electric  field  is  of  the  form  given  by 
Scharfetter  and.  Gummel:1 


v(E)  = 


(  1 

i  _ £.  +  -  - 

\  _2  k  E  +  k 
Er  34 


+  k 


(7-2) 


Values  of  the  coefficients  k^  which  are  appropriate  for  holes  in  room 
temperature  Si  were  used.  The  resulting  v(E)  has  a  low-field  mobility  of 


48o 


cm 


2  V  1  s  1  and  a  saturated  velocity  of  1.2  x  107  cm/s;  thp  complete 


function  is  plotted  in  Fig.  7.1.  The  diffusion  coefficient  was  taken  to  be 


1.  Scharfetter,  D.  L.  and  Gummel,  H.  K.,  "Large-Signal  Analysis  of  a  Silicon 
Read  Diode  Oscillator,"  IEEE  Trans,  on  Electron  Devices,  vol.  ED-l6, 

No.  1,  pp.  64-77,  January  1969"* 
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FIG.  7.1  DEPENDENCE  OF  DRIFT  VELOCITY  UPON  ELECTRIC  FIELD  FOR  ELECTRONS  AND 
HOLES  IN  Si. 
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independent  of  field  and  equal  to  12.5  cra2/s  so  that  the  Einstein 
relationship  is  satisfied  at  low  field. 

Figure  7*2  shows  profiles  of  electric  field,  hole  concentration  and 
particle  current  density  at  u>t  =  90  degrees,  when  the  total  voltage  is  a 
maximum,  for  an  operating  point  of  Jdc  =  500  A/cm2  and  =  1  V  at  12  GHz. 
For  the  remainder  of  this  section  the  term  profile  denotes  a  function  of 
space  at  a  specific  time,  e.g.,  a  doping  profile,  whereas  the  term  waveform 
denotes  a  function  of  time  only,  e.g.,  a  voltage  waveform.  From  Tig.  7.2 
it  is  apparent  that  the  depletion  width  is  approximately  5-^  4m  long.  The 
hole  concentration  is  equal  to  the  background  doping  in  the  remaining  5.6  qm 
of  the  active  region.  This  situation  is  obtained  from  the  direct  solution 
cf  the  hole  transport  equation,  not  from  the  imposition  of  any  additional 
constraint. 

Also  in  Fig.  7-2  it  is  seen  that  the  electric  field  decreases  linearly 
in  the  depletion  region  at  a  rate  determined  essentially  by  the  background 
doping.  In  the  adjacent  bulk  region  the  field,  although  it  appears  to  be 
zero  compared  to  that  in  the  depletion  region,  is  actually  uniform  at  1.3 
kv/cm.  This  is  the  value  required  to  maintain  constant  total  current  in  the 
finite  conductivity  bulk  region. 

Also  shown  in  Fig.  7*2  is  the  total  particle  current  density  which  is 

defined  as  the  sum  of  the  electron  and  hole  current  densities.  4  tie  drift 

region  the  total  particle  current  density  is  essentially  the  same  as  the 

hole  current  density  since  the  electron  current  density  is  equal  ti  J 

ns  8.x 

The  particle  current  consists  of  drift  and  diffusion  components  which  are 
determined  by  the  local  values  of  hole  and  electron  concentration  and  the 


-70- 


400 


m3  330llHVd 


D 

n 

n 


o 

o 

in 


o 

m 

«• 


i - 1 - r 

‘Q13I3  01810313 


O 

O 

ro 


O 

O 

CM 


O 

o 


o 


transport  coefficients  v(E)  and  D(E).  For  later  reference  the  induced 
current  J.  ,  is  defined  as  the  space  average  of  the  particle  current. 


1 

7 


Jp(x,t)  dx  , 


(7.3) 


and  is  a  function  of  time  only.  Because  the  scales  of  Fig.  7*2  were  chosen 
to  show  the  full  range  of  E  and  p  values  and  because  the  operating  point 
under  consideration  is  essentially  small  signal,  the  profiles  for  other 
time  instants  appear  very  similar  to  Fig.  7*2. 

Figure  7*3  shows  the  time- domain  waveforms  obtained  for  this 
operating  point.  The  RF  voltage  is  taken  to  be  a  1-V  sinusoid.  The 
resulting  current  waveforms  are  seen  to  be  approximately  sinusoidal  also. 

The  injected  current  is  defined  as  the  particle  current  at  the  interface 
of  the  avalanche  and  drift  regions,  i.e.,  x  =  1  urn  in  this  case.  The 
injected  current  is  seen  to  lag  the  voltage  by  8l  degrees.  The  induced 
current  lags  the  injected  current,  by  177  degrees.  The  device  admittance 
determined  from  the  total  current  (which  is  the  sum  of  the  induced  current 
and  the  capacitive  current  for  the  entire  active  region)  is  Y  =  -15. 5  +  jl27 
mho /cm2. 

In  Fig.  7.3  the  induced  current  was  seen  to  lag  the  injected  curient 
by  an  angle  much  larger  than  one  half  of  the  effective  transit  angle  of  the 
depletion  region.  This  is  because  the  particle  current  in  the  bulk  region 
contributes  to  the.  induced  current.  For  the  sake  of  comparison  another 

4.  X 

n  pp‘  diode  was  modeled;  it  is  identical  in  all  respects  with  the  previous 
one  except  that  the  p-layer  thickness  was  taken  to  be  5  dm  so  that  the  diode 
is  punched  through  at  breakdown  (see  Fig.  7-2).  The  computed  current 
waveforms  for  this  case  are  plotted  in  Fig.  7-^*  '"he  injected  current  is 
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FIG.  7.3  WAVEFORMS  OF  INJECTED  AND  INDUCED  CURRENT  DENSITY  FOR  A  9~um  Si  DIODE. 

(OPERATING  POINT:  =  1  V,  FREQUENCY  =  12  GHz  AND  J  =  500  A/cm2) 
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seen  to  lag  the  voltage  by  75  degrees.  The  induced  current  lags  the 
injected  current  by  60  degrees  which  i^  one  half  of  the  effective  transit 
angle  of  the  depletion  region.  The  computed  device  admittance  is  Y  = 

-55*5  +  jl26  mho/cm2.  Observe  that  the  negative  conductance  is  larger  by 
a  factor  of  two  than  for  the  unswept  epitaxial  case. 

The  difference  between  the  device  admittances  calculated  for  these 
two  cases  is  due  to  the  presence  of  ^  urn  of  bulk  p-region  in  the  former  case 
which  acts  essentially  as  a  series  resistance.  If  the  low-field  mobility  and 

a  hole  concentration  equal  to  the  background  doping  are  used,  a  resistance  of 
l.cA  cm2  is  found.  When  this  value  is  placed  in  series  with  the  device 
admittance  computed  for  the  5-pm  diode,  a  total  admittance  of  Y  =  -17  +  jl32 
is  obtained.  This  is  in  close  agreement  with  the  value  actually  computed 
for  the  9-dtn  diode. 

The  point  of  this  comparison  is  to  shjw  that  the  computer  program  can 
handle  nonpunch-through  diodes  by  extending  the  active  region  into  the  bulk. 
In  this  case  the  answer  obtain .d  shows  that  the  bulk  region  has  the  effect 
of  a  series  resistance,  in  agreement  with  physical  intuition.  In  cases  of 
large  voltage  modulation  this  approach  is  fel1.  to  be  superior  to  simply 
adding  on  an  effective  resistance. 

7»5  Silicon  Abrupt-Junction  Diodes  with  5  x  1015  Background  Doping. 

4*  +  +  + 

For  this  case  the  structures  under  consideration  are  p  nn  and  n  pp  Si 
diodes  with  a  uniform  background  doping  of  5  x  1015  cm  3  in  the  lightly  doped 
region.  Further,  the  width  of  the  lightly  doped  region  is  5  um  in  both  cases 
so  that  the  two  diode  structures  are  complementary.  Both  diodes  are  punched 
through  prior  to  breakdown. 
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To  facilitate  the  comparison  with  other  diode  analyses,  to  be 
described  later,  diffusive  transport  v?s  nf fleeted  and  the  velocities  of 
electrons  and  holes  were  taken  to  be  constant,  independent  of  the  field 
and  equal  to  107  cm/s. 

7«3»1  Small-Signal  Analysis.  An  analysis  has  recently  been 
developed  (see  Section  6)  which  calculates  the  "exact"  small-signal 
admittance  for  a  given  diode  structure.  It  is  restricted  to  small  signal 
since  it  performs  linear  perturbations  upon  the  static  current  and  field 
profiles;  these  are  obtained  from  a  static  analysis  described  previously.2 
However,  it  requires  no  a  priori  assumption  about  the  localization  of  the 
impact  ionization  and  it  handles  unequal  ionization  rates  and  drift  velocities 
for  electrons  and  holes  (the  drift  velocities  are  assumed  constant). 

Figure  7-5  shows  the  small-signal  admittance  curves  obtained  from 
the  exact  analysis  for  the  complementary  Si  diodes  biased  at  50C,  1000  and 
2000  A/cm2.  It  is  seen  that  at  each  bias  current  level  the  peak  negative 
conductance  is  about  50  percept  larger  for  the  p-type  han  for  the  n-type 

V 

diode.  This  is  due  to  the  inequality  of  the  ionization  rates,  which  results 
in  a  narrower  avalanche  region  for  the  p-type  structure  as  has  been  discussed 
previously.3 

In  order  to  model  the  complementary  Si  diodes  with  the  large-signal 
analysis  of  this  section,  it  is  necessary  to  select  the  effective  avalanche 


2.  Schrocder,  W.  E.  and  Haddad,  G.  I.,  "Avalanche  Region  Width  in  Various 
Structures  of  IMPATT  Diodes,"  Proc.  IEEE  (Correspondence),  vol.  59> 

No.  8,  pp.  1245-1248,  August  1971. 

3-  Gilden,  M.  and  Hines,  M.  E.,  "Electronic  Tuning  Effects  in  Read 
Microwave  Avalanche  Diodes, "  IEEE  Trans .  on  Electron  Devices, 
vol.  ED-13,  No.  1,  pp.  169-175,  January  19667 
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FIG.  7-5  SMALL  -SIGNAL  ADMITTANCE  DETERMINED  FROM  THE  EXACT  ANALYSIS  FOR  COMPLEMENTARY  Si  DIODES. 
(BACKGROUND  DOPING  =  5  x  1015  cm"3  AND  WIDTH  =  5  |im) 


region  width  appropriately.  Previously,  using  a  somewhat  arbitrary 
definition  of  avalanche  width,  i.e.,  the  region  where  95  percent  cf  the 
particle  current  is  generated,  effective  avalanche  widths  of  IT  and  ^0 
percent  were  established  for  p-  and  n-tvpe  diodes  of  this  doping  level, 
respectively.  In  view  of  the  arbitrariness  of  the  95  percent  criterion, 
it  was  decided  to  make  a  direct  comparison  of  the  small -signal  admittance 
calculated  from  the  exact  analyois  with  that  obtained  from  the  localized 
avalanche  model  using  different  effective  avalanche  widths.  Since  for  the 
case  of  constant  velocity  and  no  diffusion  the  localized  avalanche  model 
reduces  to  the  familiar  Read  model,  an  analytic  expression  is  obtained  for  the 
small-signal  admittance:3 


by  the  avalanche  width  as  well  as  tne  ionization  rates. 

4*  4* 

Figure  7-6  shows  this  comparison  for  the  p  nn  diode  biased  at 
500  A/cm2.  Observe*  that  the  exact  admittance  curve  does  not  pass  through 
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1|  FIG.  7.6  COMPARISON  OF  EXACT  SMALL-SIGNAL  ADMITTANCE  WITH  THE  APPROXIMATE 

i  DIODE  EXPRESSION  FOR  A  5-nm  p+nn+  Si  DIODE  BIASED  AT  500  A/cm2. 
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the  origin  of  the  admittance  plane.  The  conductance  and  susceptance  change 
sign  at  approximately  6.6  GHz  and  7*1  GHz,  respectively.  On  the  other  hand, 
the  Read  model  admittance  curves  all  pass  thro'igh  the  origin  at  f  =  f  ,  as 
may  toe  seen  from  inspection  of  Eq.  7;^«  Thus,  no  matter  what  choice  of 
t^/t  and  f  is  made,  the  curves  cannot  toe  in  complete  agreement  near  the 
avalanche  frequency.  However,  the  frequencies  which  are  generally  of  the 
most  interest  are  those  where  peak  negative  conductance  occurs;  here  the 
agreement  can  toe  made  excellent  as  shown  in  the  figure. 

In  Fig.  7*6  the  avalanche  frequency  is  taken  to  be  7*1  GHz  for  all 
of  the  Read  model  curves.  It  is  seen  that  the  effect  of  widening  the 
avalanche  region  is  to  reduce  the  peak  negative  conductance  substantially 
and  to  increase  somewhat  the  frequency  for  peak-G  concomitant  with  the 
reduced  drift  region  length.  Although  it  is  not  shown  in  the  figure,  the 
effect  of  increasing  avalanche  frequency  while  maintaining  fractional 
avalanche  width  constant  is  to  increase  the  peak  negative  conductance.  Of 
course,  the  frequency  at  which  the  admittance  zero  occurs  is  equal  to  the 
avalanche  frequency. 

It  is  seen  in  Fig.  7»6  that  a  choice  of  l^/l  =  0.75  gives  good 
agreement  between  the  Read  admittance  curve  and  the  exact  admittance, 
particularly  near  the  peak  negative  conductance  frequencies.  This  corresponds 
to  a  fractional  avalanche  width  of  25  percent  or  1.25  pm.  A  similar  comparison 
was  made  for  this  diode  operated  at  2000  A /cm2  and  again  good  agreement  was 
obtained  with  the  exact  curve  using  an  avalanche  frequency  increased  by  a 
factor  of  two  (square  root  of  the  current  ratio)  and  an  avalanche  width  of 
25  percent.  Thus  x  =  1.25  pm  is  used  in  the  remainder  of  this  section  for 
the  p+nn+,  5  x  1015  cm  3  Si  diode.  Similar  comparisons  were  carried  out  for 
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the  complementary  n  pp  diode  and  yielded  a  narrower  effective  avalanche 
region  equal  to  15  percent  of  the  total  depletion  width  or  0.75  pm;  this 
value  is  used  subsequently  for  the  large-signal  calculations. 

7.J>.2  Large -Signal  Analysis.  Recently,  a  large-signal  diode 
analysis  lias  been  developed  in  connection  with  a  study  of  the  high-efficiency 
mode,  which  also  may  be  used  for  studying  the  IMPA1T  mode.4  It  permits 
distributed  ionization  and  requires  saturated  velocities  as  does  the  exact 
small-signal  analysis.  Under  snail  drive  conditions  it  computes  admittance 
values  in  good  agreement  with  those  from  the  exact  small-signal  analysis; 
since  the  programs  are  independent  this  agreement  serves  as  a  validity  check. 
It  is  of  interest  to  compare  the  results  of  this  program  with  those  from  the 
localized  avalanche  model  of  this  section.  For  ease  of  identification  the 
models  will  be  referred  to  as  distributed  and  Read,  respectively. 

For  the  remainder  of  Section  7*5«2  the  operating  point  of  thu 
complementary  Si  diodes  will  be  restricted  to  a  bias  current  of  1200  A/cm2 
and  an  operating  frequency  of  1^  GHz,  which  is  near  optimum  at  large  signal. 
Figure  7«7  shows  the  device  admittance  as  a  function  of  the  RF  voltage  level 
computed  for  the  Read  model  for  the  complementary  Si  diodes.  It  is  seen 
that  at  small  signal  the  negative  conductance  for  the  n+pp+  diode  ?s  about 
50  percent  larger  than  for  the  complement,  that  the  negative  conductance 
decreases  monetonically  with  voltage  for  both  and  that  at  50  V  the  p-type 
has  a  negative  conductance  only  20  percent  larger  than  the  n-type.  The 
susceptance  for  both  diodes  increases  with  voltage,  approaching  the  depletion 
layer  capacitance  value  at  large  drive  levels. 

4.  Haddad,  G.  I.  et  al.,  "Microwave  Solid-State  Device  and  Circuit  Studies," 
Quarterly  Progress  Report  No.  5,  Contract  No.  F30602-71-C-0099>  Electron 
Physics  Laboratory,  The  University  of  Michigan,  Ann  Arbor,  June  1971* 
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FIG.  7-7  DEVICE  ADMITTANCE  AS  A  FUNCTION  OF  RF  LEVEL  FOR  THE  COMPLEMENTARY 

Si  DIODES.  (BACKGROUND  DOPING  =  5  x  1015  cm"3  AND  WIDTH  =  5  4m. 

OPERATING  POINT:  FREQUENCY  =  lb  GHz  AND  J  =  1200  A /cm2) 
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Figure  7 *8  compares  the  induced  current  waveforms  computed  from 
the  distributed  and  Read  models  for  the  n+pp+  diode  at  V™,  =  50  V-  They 

IYT 

Eire  seen  to  be  in  very  good  agreement  and  consequently  the  diode  admittance 
valuer,  are  also.  The  Read  model  admittance  is  Y  =  -15-9  +  jl7^  and  the 
distributed  model  value  is  Y  =  -l6.3  +  jl77*6  mho/cm2.  In  Fig.  7»®  it  is 
seen  that  the  agreement  is  very  good  for  the  portion  of  time  when  the 
induced  current  rises  from  a  low  value  to  the  peak  value;  this,  of  course, 
is  the  time  when  avalanche  generation  of  carriers  takes  place.  The  hump 
in  induced  current  slightly  before  l8o  degrees  is  due  to  the  contribution 
of  electron  current.  After  about  200  degrees  the  induced  current  is  constant 
because  all  the  electrons  have  been  collected  and  the  hole  pulse  is  drifting 
at  constant  velocity.  Slightly  before  360  degrees  the  hole  pulse  begins  to 
be  collected  at  the  far  end  of  the  drift  region. 

Figure  7*9  shows  tne  particle  current  profile  at  several  times  in  the 
RF  cycle  calculated  from  the  Read  model  for  the  same  operating  point  as 
Fig.  7.8.  The  four  time  instants  shown  correspond  to  the  generation  and 
launch  of  the  hole  current,  pulse.  Figure  7-10  compares  the  particle  current 
profiles  at  cot  =  l8o  degrees  obtained  from  the  Read  and  distributed  models. 

In  the  dis  ributed  model  the  particle  current  is  found  to  be  nonuniform  in 
the  avalanche  .region.  Also,  the  hole  pulse  in  the  distributed  model  lags 
that  of  the  Read  model  by  0.25  urn,  which  corresponds  to  the  hole  pulse  being 
launched  from  the  middle  of  the  effective  avalanche  region. 

Figure  7*11  compares  the  induced  current  waveforms  calculated  from 
the  two  models  for  the  p  nn  diode  at  V^p  =  ;0  V,  i.e.,  the  case  complementary 
to  Fig.  7.8.  The  waveforms,  are  seen  to  be  very  similar  in  shape  but  the 
induced  current  from  the  distributed  model  lags  the  other  by  about  15  degrees. 
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FIG.  7.8  COMPARISON  OF  INDUCED  CURRENT  WAVEFORMS  DETERMINED  BY  THE  LARGE- 
SIGNAL  MODELS  FOR  A  5-Uir.  n+pp+,  5  x  1015  cm'3  DIODE.  (OPERATING 
POINT:  Vgj,  =  50  V,  FREQUENCY  =  14  GHz  AND  Jdc  =  1200  A /cm2) 
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7*9  PARTICLE  CURRENT  PROFILES  AT  SEVERAL  TIMES  IN  THE  RF  CYCLE  FOR  A 

5-iJm  n  pp  ,  5  x  1015  cm”3  DIODE.  (OPERATING  POINT:  VnT?  =  50  V, 

Kb 

FREQUENCY  =  1^  GHz  AND  .TJc  =  1200  A/cm2) 
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FIG.  7.10  COMPARISON  OF  PARTICLE.  CURRENT  PROFILES  AT  cJ.  =  3 8o  DEGREES 

DETERMINED  BY  THE  LARGE-SIGNAL  MODELS  FOR  A  b-pai  n+pp+,  5  x 

cm"3  DIODE.  (OPERATING  POULT;  V  =  V,  FREQUENCY  =  'J‘  C-E 

AND  J  ,  =  1200  A /cm2) 
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As  a  result  the  admittances  do  not  agree  as  closely  as  in  the  complementary 
case.  The  admittances  computed  are  Y  =  -15.^  +  jl68  mho/cm2  and 
Y  =  -17*9  +  J171  mho/cm2  for  the  Read  and  distributed  models,  respectively. 
Again,  the  hump  in  induced  current  near  l8o  degrees  is  due  to  "minority" 
earlier  current,  in  this  case  holes.  The  hump  is  more  pronounced  than  in 
Fig.  7*8  because  the  avalanche  region  is  wider. 

As  we  have  already  seen  in  the  small-signal  case,  negative  conductance 

decreases  with  increasing  avalanche  width.  In  the  two  large -signal  cases 

just  described  smaller  negative  conductance  was  calculated  from  the  Read 

model  for  25  percent  than  for  15  percent  effective  avalanche  widths.  This 

appears  to  be  a  general  property  of  the  Read  model  as  is  shown  in  Fig.  7*12; 

here  large-signal  admittances  are  shown  which  were  calculated  for  a  diode 

of  5  uni  total  length  and  various  avalanche  widths.  The  operating  point  for 

all  cases  is  =  50  V,  frequency  =  l4  GHz  and  J,  =  1200  A/cm2.  The 
Kr  a'. 

avalanche  widths  other  than  15  and  25  percent  do  not  necessarily  correspond 
to  any  particular  diode  structure  but  rather  are  used  to  show  the  effect  of 
the  avalanche  width  parameter  in  the  model. 

Table  7«1  summarizes  the  result:,  obtained  from  the  two  models  for 
the  complementary  Si  diodes  at  the  1^  GHz,  50  V  RF,  1200  A/cm2  operating 
point.  For  the  p-type  diode  the  two  models  give  results  which  are  in 
excellent  agreement;  the  calculated  efficiency  is  slightly  less  than 
20  percent.  For  the  n-type  diodes  the  two  models  differ  somewhat  as 
described  previously.  However,  what  may  be  more  important  than  the 
difference  is  that  either  model  predicts  that  the  efficiencv  is  not  greatly 
different  from  that  of  the  complementary  diode,  e.g.,  two  percent  less  in 
the  Read  model  vs.  two  percent  more  in  the  distributed  model.  Two 
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qualifications  should  he  noten  with  respect  to  the  efficiency  values  quoted 
here.  First,  the  operating  point  is  not  necessarily  that  which  yields  peak 
efficiency  for  either  structure.  Second,  the  efficiencies  quoted  are 
optimistic  since  in  both  models  the  field  in  the  drift  region  was  computed 
to  be  very  small  near  cot  =  270  degrees  which  would  lead  to  a  smaller  induced 
current  if  the  velocities  were  permitted  to  become  unsaturated. 


Table  7-1 

Large-Signal  Data  for  Complementary  Diodes 


Type 

Model 

G  (mho /cm2) 

B(mho/cm2) 

vac  <v> 

PRj/kW/cm2) 

1  W 

+  + 
n  pp 

Read 

-15.9 

174 

84.5 

19.9 

19.7 

+  + 
n  pp 

distributed 

-16.3 

177.6 

85.4 

20.4 

19.9 

+  + 
p  nn 

Read 

-13.4 

167.5 

76.5 

16.7 

18.2 

+  + 
p  nn 

distributed 

-17.9 

170.8 

• 

QO 

22.4 

22.3 

7.4  Millimeter -Wave  Silicon  Diodes .  The  diode  structure  under 

■f  + 

consideration  in  this  section  is  a  one-sided,  p  nn  Si  junction  with 
n-region  doping  of  7  x  1016  cm  3  and  width  of  0.50  pm.  The  dc  current  is 
taken  equal  to  104  A/cm2  in  all  cases.  Material  parameters  appropriate  to 
Si  at  200 °C  are  used  throughout. 

The  small-signal  admittance  of  this  structure  was  evaluated  using 
the  exact  analysis  mentioned  in  the  previous  section.  For  the  purpose  of 
this  analysis  the  velocities  of  electrons  and  holes  were  taken  to  be 
constant  and  equal  to  8.8  x  10s  cm/s.  Since  the  diode  is  punched  through 
at  breakdown  with  a  minimum  field  of  9*8  x  104  v/cm,  this  is  a  reasonable 
approximation  for  the  static  and  small-signal  case.  The  ionization  rates 
were  taken  to  be  of  the  form  A[exp(-b/E)  1,  where  for  electrons 
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A  =  1.08  x  10s  cm  1  and  b  =  1.34  x  106  v/cm  and  for  holes  A  =  2.56  x  106  cm 
and  b  =  2. 44  x  10s  v/cm.  These  values5  are  bei, ieved  to  be  appropriate  to 
Si  at  high  fields  (~  6  x  ]05  v/cra)  and  200°C.  Figure  7.13  shows  the  ex> ct 
small-sigral  admittance  which  was  computed  using  these  parameters.  The 
frequencies  for  zero  conductance  and  susceptance  are  35  and  42  GHz, 
respectively.  The  breakdown  voltage  was  found  to  be  16.9  V  with  low 
current  ana  If. 2  V  at-  the  operating  current  of  104  A/cm2. 

The  analytic  small-signal  admittance  expression  of  Eq.  7*4  was 
evaluated  for  several  choices  of  avalancne  width  and  avalanche  resonant 
frequency.  A  relative  avalanche  width  of  30  percent  (or  0.150  pm)  and  an 
avalanche  frequency  of  38  GHz  were  found  to  give  a  good  fit;  the  resulting 
admittance  is  also  plotted  in  Fig.  7.13*  Observe  that  the  effective 
avalanche  width  is  wider  than  for  the  n-type  X-band  Si  diode;  this  is 
because  the  ionization  rates  are  flatter  at  high  field  which  leads  to  a 
spread  out  generation  region.3 

Using  the  large-signal,  localized  avalancne  computer  program,  values 
for  device  admittance  were  computed  as  a  function  of  frequency  and  voltage 
level  and  are  shown  in  Fig.  7*l4.  For  this  program  the  drift  velocity  of 
electrons  was  taken  to  have  the  field  dependence  sho' n  in  Fig.  7.1;  the 
saturated  value  is  8.8  x  10s  cm/s  and  the  low-field  mobility  is  220  cm2/v-s. 
The  diffusion  coefficient  was  taken  to  be  field  independent  and  equal  to 
9  cm2/s  in  agreement  with  the  Einstein  relationship  at  low  field.  An 
effective  ionization  rate  was  determined  which  gives  the  desired  breakdown 
voltage  and  avalanche  resonant  frequency;  it  is  close  to  that  of  electrons. 


5.  Bowman,  L.,  private  communication. 
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In  Fig.  7.2>  the  VDT,  =  0.1^  curve  corresponds  to  small  signal  since 
the  computed  ac  current  was  found  to  he  sinusoidal  and  very  small  compared 
to  the  dc.  It  is  in  good  agreement  with  the  small-signal  admittance  shewn 
in  Fig.  7.1?  which  was  computed  from  Eq.  7«^,  except  that  the  neaiv  negative 
conductance  here  is  slightly  smaller.  This  difference  is  attributed  to  the 
effects  of  unsaturated  velocity  and  diffusion  which  were  not  incorporated 
in  Eq.  7.^-.  The  curve  for  V™  =  1.^  V  is  seen  to  be  near  the  beginning  of 

Kr 

the  nonlinear  regime  because  it  differs  only  slightly  from  that  for 

V™  =  0.l4  v.  For  voltages  above  7*1  V  data  were  taken  only  for  f  =  8o  GHz. 

Kr 

Figure  7*15  shows  the  power  and  efficiency  calculat  >d  at  f  =  8o  GHz 
as  a  function  of  RF  voltage.  It  is  seen  that  power  and  efficiency  saturate 
at  approximately  10  V  and  then  decrease  with  increased  drive  level.  The 
peak  efficiency  calculated  is  l4,l  percent.  The  power  and  efficiency 
saturate  because  the  induced  current  waveform  develops  a  dip  at  approximately 
270  degrees.  The  size  of  the  dip  increases  with  voltage  level  as  shown  in 
Fig.  7.18.  Recause  the  dip  occurs  at  the  time  when  the  RF  voltage  is 
minimum  it  seriously  degrades  the  generated  RF  power. 

The  cause  of  the  dip  is  shown  in  Fig.  7.17  which  gives  the  electric 
field  and  electron  concentration  profiles  for  several  time  instanbs  at  the 
l4-v  operating  point.  It  is  seen  that  the  electric  field  in  the  drift 
region  becomes  negative  which  results  in  the  electrons  drifting  opposite 
to  the  normal  direction.  Also,  electrons  flow  into  the  drift  region  from 
the  contact  which  normally  collects  electrons,  in  what  is  essentially 
depletion  lajer  width  modulation.  Both  of  these  processes  give  negative 
particle  current  which  is  manifested  as  the  dip  in  induced  current. 
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FIG.  7-15  POWER  AND  EFFICIENCY  AS  A  FUNCTION  OF  RF  LEVEL  FOR  A  0.5- urn 
p+nn+  Si  DIODE.  (OPERATING  POULT:  FREQUENCY  =  80  GHz  AND 
J.  =  104  A/cm2) 
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7*5  Conclusions  and  Program  ior  the  Next  Quarter.  In  the  preceding 
sections  some  results  which  have  been  obtained  from  the  large-signal 
localized  avalanche  computer  program  have  been  presented  and  discussed.  It 
has  bean  shown  that  the  program  can  treat  diode  structures  in  which  a 
portion  of  the  active  region  is  undepleted.  A  method  has  been  described 
for  establishing  an  effective  avalanche  width  by  fitting  an  exact  small- 
signal  admittance  curve.  Large-signal  results  for  an  X-band  Si  diode  using 
this  effective  avalanche  width  were  found  be  in  good  agreement  with 
another  large-signal  program  which  permits  distributed  avalanche.  Results 
have  also  been  obtained  for  a  millimeter-wave  Si  diode  which  show  that  power 
and  efficiency  saturate  with  voltage  drive  due  to  physical  processes  which 
occur  in  the  drift  region. 

During  the  next  period  this  analysis  will  be  used  to  study  other 
diode  structures  and  in  particular  to  study  GaAs  diodes.  The  effect  of 
lattice  temperature  upon  device  performance  will  also  be  studied  in  a 
systematic  way. 

8.  Solid-State  Device  Fabrication 
Supervisor:  N.  A.  Masnari 
Staff:  J.  East, 

8.1  Introduction.  The  objective  of  this  phase  of  the  program  is  to 
fabricate  IMPATT  diodes  using  the  facilities  of  the  Electron  Physics 
Laboratory.  Significant  progress  has  been  made  toward  this  goal  during  this 
reporting  period  with  devices  being  fabricated  from  three  different  types 
of  wafers : 

1.  nn+  epitaxial  material  in  which  the  p+  layer  is  formed  by 
diffusion. 
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2.  Double  epitaxial  p+n  layers  on  n+  substrate. 

3.  pp+  epitaxial  material  with  the  n+  layer  formed  by  ion  implantation 
of  phosphorus  ions. 

The  introduction  of  integral  gold-plated  heat  sinks  has  resulted  in 
significant  improverent  in  the  quality  of  the  packaged  diodes.  A  large 
number  of  diodes  have  exhibited  good  thermal  resistances  as  well  as  powers 
in  excess  of  600  mW  and  efficiencies  as  high  as  6  percent. 

3.2  Diode  Fabrication.  Changes  in  the  fabrication  process  have  been 
made  to  improve  the  thermal  resistance  of  the  diodes.  Using  the  process 
described  in  Reference  1,  the  wafer  is  cleaned,  a  junction  is  produced  by 
diffusion  or  ion  implantation  and  the  wafer  is  thinned  and  metallized 
(chrome-gold).  It  is  at  this  point  that  significant  improvement  has  been 
achieved  by  gcld-plating  the  wafer  in  a  commerical  plating  solution  to 
produce  an  integral  heat  sink.  The  thickness  of  the  gold  placing  has  varied 
from  0.001  inch  to  0.003  inch. 

Several  - oehniques  have  been  used  in  the  plating  process.  One  method 
is  to  plate  the  gold  over  the  entire  junction  side  of  the  wafer.  The 
substrate  side  is  then  covered  with  an  array  of  small  diameter  photoresist 
circles  after  which  the  metal  not  under  the  dots  is  etched  away.  The  wafer 
is  then  submitted  to  a  silicon  etch  which  removes  all  of  the  Si  except  that 
under  the  dots  so  that  the  final  result  is  an  array  of  Si  mesas  attached 
to  the  gold-plated  layer.  The  excess  gold  rims  projecting  over  the  Si  mesas 
(a  result  of  undercutting  during  etching)  are  removed  by  subjecting  the  wafer 


1.  East,  J.  and  Masnari,  N.  A.,  "IMPATT  Diode  Fabrication,"  EPL  Memo 
No.  71-l-°37250,  Electron  Physics  Laboratory,  The  University  of 
Michigan,  Ann  Arbor,  2  June  1971. 
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to  the  gpld  etch  once  again  Since  the  tops  of  the  gold  caps  are  protected 
by  the  photoresist,  the  effect  is  for  the  etch  to  attack  only  the  gold  rims 
from  underneath.  This  eliminate.,  the  troblems  experienced  in  the  past  of 
the  gpld  rims  folding  over  and  shorting  the  substrate  and  .-function  sides  of 
the  diode.  The  wafer  is  then  separated  into  individual  diodes  by  cutting  the 
gold  plating  with  a  sharp  knife  or  razor  blade. 

A  second  method  is  to  use  photolithographic  techniques  to  form  a 
pattern  of  0.015-inch  diameter  windows  in  the  layer  of  photoresist  on  the 
Junction  side  of  the  wafer.  The  wafer  is  then  waxed  to  a  piece  of  glass  to 
protect  the  back  and  to  give  mechanical  support.  When  the  wafer  is  inserted 
into  the  plating  solution  only  the  open  windows  are  plated.  After  being 
plated  uhe  diodes  are  formed  on  the  substrate  side  as  described  above  with 
the  additional  precaution  that  the  diodes  must  be  aligned  with  the  gold- 
plated  circles  on  the  junction  side.  In  this  case  +l,e  diodes  are  separated 
from  each  other  after  the  final  silicon  etch.  The  individual  diodes  are 
then  cleaned  in  an  ultrasonic  bath  sequentially  in  trichlorethylene,  acetone 
and  Freon  after  which  they  are  baked  at  150°C  overnight. 

The  finished  diodes  vary  in  diameter  from  0.005  inch  to  0.005  inch. 

The  first  process  described  above  results  in  the  diodes  being  mounted  on 
square  gold-plated  pads.  In  the  second  process  the  diodes  are  located  on 
0.015-incn  diameter  gold-plated  cylinders.  Each  unit  is  then  thermo¬ 
compression  bonded  into  a  package  by  bringing  the  bonding  tip  down  around 
the  periphery  of  the  diode,  being  careful  not  to  touch  the  diode  itself. 

A  0.005-inch  diameter  tip  is  used  with  75,000  psi.  The  heavy  pressure 
ensures  a  good  thermal  contact  between  the  package  rnd  the  integral  gold- 
plated  heat  sink.  There  is  no  danger  of  damaging  the  diode  at  this  point 
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because  it  is  not  touched  by  the  bonding  tip.  After  the  gold  is  bonded,  a 
much  lighter  pressure  is  used  to  bond  a  wire  from  the  top  of  the  diode  to 
the  upper  flange  of  the  package. 

8.3  Experimental  Results .  A  group  of  3^  ion-implanted  diodes  was 
fabricated  using  the  above  techniques.  The  thermal  resistances  of  the 
diodes  with  good  breakdown  characteristics  were  measured  using  the  standard 
technique  and  the  definition  of  thermal  resistance  as  described  by 
Haddad  et  al.2  The  average  thermal  resistance  was  21.7°c/w  with  several 
being  as  low  as  17°c/w.  In  some  cases  the  diodes  were  able  to  withstand 
13.5  W  of  dc  power  before  burnout  occurred. 

A  second  group  of  19  diodes  was  fabricated  using  a  double  epitaxial 
wafer.  This  group  had  an  average  thermal  resistance  of  2^°c/w.  Part  of 
this  higher  thermal  resistance  can  be  explained  by  the  fact  that  the 
junction  was  much  deeper  in  the  material. 

From  the  above  results  it  appeared  that  the  thermal  resistance  had 
become  a  secondary  problem  so  three  more  wafers  were  fabricated  in  the  hope 
of  improving  the  RF  operation  of  the  diodes.  Two  of  the  wafers  were  from 
single  epitaxial  layer  material  with  a  3.5-pm  thick,  0.6  fi-cm  epitaxial 
layer.  One  of  these  (wafer  A)  had  a  0.5-pm  deep  dif^asion  while  the  other 
(wafer  B)  had  a  0.7-mei  deep  diffusion.  The  third  wafer  (wafer  C)  was  a 
double  epitaxial  wafer  (p+nn+)  with  the  n -layer  having  a  carrier  concentration 
of  approximately  1016  cm  3.  The  breakdown  voltage  of  these  diodes  is 
approximately  60  to  65  V.  Although  all  of  the  diodes  have  not  bsen  packaged 
and  tested,  several  of  them  from  each  wafer  have  been  tested  for  RF 

2.  Haddad,  G.  I.,  Greiling,  P.  T.  and  Schroeder,  W.  E.,  "Basic  Principles 
and  Properties  of  Avalanche  Transit-Time  Devices, "  IEEE  Hans,  on 
Microwave  Theory  and  Techniques,  vol.  MTl-l8,  No.  11,  pp.  752-772, 

November  1970 
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performance.  Diodes  from  wafer  A  have  averaged  approximately  100  mW  CW 
power  output  in  the  6.5  to  8.5  GHz  frequency  range.  Initial  tests  of  diodes 
from  wafer  B  averaged  above  150  mW  in  the  same  frequency  range.  In  both 
cases  the  bias  currents  were  approximately  90  mA.  Recently,  two  diodes 
from  wafer  B  have  yielded  in  excess  of  600  mW  at  6.9  GHz.  However,  these 

diodes  were  rather  noisy  at  these  low  frequencies  while  at  higher  frequencies 
they  were  unable  to  generate  more  than  a  few  hundred  milliwatts. 

Diodes  from  wafer  C  have  yielded  power  outputs  ranging  from  100  to 
360  mW  in  the  7-7  to  11  GHz  frequency  range.  A  listing  of  some  of  these 
diodes  is  shown  in  Table  8.1.  It  should  be  mentioned  that  several  of  these 
diodes  burned  out  while  tuning  because  of  circuit  problems,  in  particular, 
most  of  the  difficulty  was  associated  with  movement  of  the  center  bias 
conductor  during  tuning. 


Table  8.1 

Diodes  from  Wafer  C  (Double  Ep^'+axial  Material) 

Diode  Humber  Bias  Current  (mA)  Power  Output  (mW)  Efficiency  (<&) 


1 

70 

100 

1.9 

2 

55 

205 

5.0 

3 

4-0 

110 

3.8 

k 

4o 

170 

5-7 

5 

40 

165 

5-5 

6 

50 

210 

5.6 

7 

30 

105 

4.0 

8 

50 

360 

9.8 
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Following  oreliminary  evaluation  of  the  afore-mentioned  diodes,  another 
epitaxial  wafer  (D)  having  the  same  specifications  as  A  and  B  was  diffused  to 
create  a  deep  junction.  The  diodes  fabricated  from  this  wafer  are  somewhat 
smaller  (~  0.003-inch  diameter)  than  the  others.  At  the  same  time  as  these 
diodes  were  being  fabricated,  an  improved  copper  collet  and  heat  sink  was  also 
constructed.  The  combination  of  new  diodes  and  improved  heat  sink  resulted  in 
a  vast  improvement  in  the  performance  of  the  diodes.  Thus  far  four  diodes  have 
been  carefully  evaluated  using  the  new  arrangement.  Each  diode  has  generated 
in  excess  of  300  mW  of  power  output  across  the  frequency  band  from  7.4  to  10.6 
GHz  for  bias  currents  of  130  mA.  Figure  8.1  illustrates  the  power  output  and 
efficiency  obtained  -from  one  of  the  diodes  as  a  function  of  frequency  for 
various  bias  currents.  The  power  output  at  each  bias  level  extends  from  7-8  to 
11  GHz.  The  maximum  power  of  480  mW  at  9-5  GHz  is  achieved  at  a  bias  current 
of  I30  mA.  The  maximum  efficiency  of  4  percent  also  occurs  at  this  point. 

Figure  8.2  illustrates  an  even  better  performance  by  a  second  diode. 
Again  the  diode  is  tunable  over  a  broad  frequency  range.  In  addition,  it  gen¬ 
erates  a  maximum  power  output  of  64o  mW  at  9-4  GHz  for  bias  current  of  130 
mA.  The  maximum  ®ff  i  H  ency  also  occurs  at  this  point  and  is  approximately  5-4 
percent.  It  is  apparent  that  both  the  power  output  and  efficiency  have  their 
greatest  incremental  increase  in  3oing  from  120  to  I30  mA  cf  bias  current. 

The  other  diodes  had  characteristic.  milar  to  the  above.  None  of 
them  have  been  pushed  far  enough  to  cause  burnout.  Based  on  the  measured 
thermal  resistances,  it  is  anticipated  that  an  additional  20  to  30  mA  of  bias 
current  should  be  possible  without  destroying  the  diodes.  Extension  of  the 
measurements  to  the  higher  bias  currents  will  be  implemented  as  soon  as  the 
fabrication  of  the  newly  designed  copper  heat  radiating  structure  is  completed. 


DIODE  1  FABRICATED  FROM  WAFER  D 


OUTPUT  POWER - EFFICIENCY 
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.Tv  ELS.  DIOLE  2  FABRICATED  FROM  WAFER  D. 


8.4  Conclusions .  The  fabrication  of  IMPATT  diodes  has  improved 


significantly  during  this  reporting  period.  In  particular,  the  incorporation 
of  the  integral  gold-plated  heat  sink  has  resulted  in  improved  thermal 
resistance.  In  addition,  the  recent  acquisition  of  improved  material  has 
also  aided  in  the  attainment  of  better  power  outputs  and  efficiencies. 

Thermal  resistances  as  low  as  17°C/W  along  with  maximum  power  outputs  of 
600  mW  and  peak  efficiency  of  6  percent  indicate  that  the  fabrication 
technique  is  coming  under  control. 

8.3  Program  for  the  Next  Quarter.  During  the  next  quarter  the 
remaining  diodes  from  wafers  A,  B,  C  and  D  will  be  packaged  and  tested.  Also, 
other  wafers  .fill  be  fabricated  to  explore  th'  relationship  between  doping 
profile  and  density  and  the  KF  operation  of  the  diodes.  Now  that  the  power 
and  efficiency  appear  reasonabl  more  effort  will  be  devoted  to  increasing 
tne  operating  frequency.  The  primary  concern  will  be  with  optimization  of 
the  circuit  to  improve  the  RF  performance  of  the  diodes.  At  the  present 
time  the  major  limitations  appear  to  be  associated  with  the  circuit  itself. 

£.  Instabilities  in  Germanium 
Supervisor:  J.  E.  Rowe 
Staff:  E.  H.  Sigmar. 

9.1  Introduction.  The  previous  quarterly  progress  report  indicated 
accuracy  problems  that  were  present  in  the  computer  subroutine  used  to 
calculate  collision  integrals  for  the  dc  model  of  germanium.  Thit 
subroutine  used  the  value  of  the  integrand  at  a  few  fixed  points  in  the 
interval  of  integration  to  evaluate  the  integral  by  Simpson's  rule. 

However,  in  some  cases  it  was  found  that  the  integrand  had  high  amplitude 
spikes  within  the  range  of  integration  which  could  cause  the  Simpson's  rule 
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integration  to  be  inaccurate  when  only  a  few  fixed  points  were  used. 
Modifications  have  been  made  to  the  computer  subroutine  to  alleviate  this 
problem. 

9.2  Modifications  in  the  Simpson ’ s  Integration  Procedure.  In  order  to 
separate  the  numerical  integration  problems  from  the  collision  integral 
subrou+ino  «  separate  multiple-entry  subroutine  has  been  written  for  the 
Simpson's  integration.  The  collision  integral  subroutine  has  been  modified 
to  call  on  this  subroutine  as  needed.  The  Simpson's  integration  subroutine 
uses  a  variable  step  size  to  achieve  the  desired  accuracy.  Since  increasing 
the  number  of  sample  points  along  the  entire  integration  interval  would  be 
an  inefficient  method  of  improving  the  accuracy,  the  subroutine  increases 
the  number  of  sample  points  only  for  the  portions  of  the  integration  interval 
where  it  is  necessary.  This  is  achieved  as  follows. 

The  subrc  tine  starts  with  five  equally  spaced  points  along  the 
integration  interval.  The  value  of  the  integral  is  calculated  from  Simpson's 
rule  for  five  points.  The  two  end  points  and  the  midpoint  of  the  interval 
are  then  used  and  the  integral  value  is  calculated  from  Simpson's  rule  for 
three  points.  The  three  and  five  point  values  are  compared  and  if  they 
are  within  the  specified  maximum  allowable  discrepancy  the  five  point  value 
is  used  for  the  integral.  If  the  three  and  five  point  values  are  not  within 
the  allowable  limits  ther  the  integrand  value  is  calculated  at  four  new 
points --the  points  midway  oetween  each  pair  of  previously  adjacent  points. 

The  integration  interval  may  now  be  divided  in  half  resulting  in  five 
equally  spaced  points  in  each  half- interval  at  which  the  integrand  value 
has  a  calculated.  The  three  and  five  point  integration  values  may  now 
be  compared  in  each  half -interval.  If  the  desired  accuracy  has  been 
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attained  for  a  half  interval,  the  five  point  value  can  be  used  for  the 
half -interval.  If  the  desired  accuracy  has  not  been  attained,  the  half¬ 
interval  can  be  further  subdivided  and  the  same  procedure  used  on  each 
subinterval.  In  this  way  extra  sample  points  are  taken  only  in  the  local 
vicinity  of  troublesome  points  rather  than  across  the  entire  integration 
interval. 

9»3  Results .  The  computer  subroutine  to  calculate  the  integral  by 
Simpson's  rule  has  been  implemented  and  the  results  have  been  found  to  be 
highly  accurate.  However  in  some  cases  it  has  been  noted  that  the  Simpson's 
integration  subroutine  uses  over  a  thousand  sample  points  to  calculate  the 
integral.  Because  of  the  numerous  integrations  performed,  this  quickly 
leads  to  excessive  computer  time  for  the  entire  program. 

One  approach  which  has  been  used  to  alleviate  the  problem  of 
excessive  computer  time  has  been  to  change  the  error  specification  for  the 
Simpson's  integration.  The  accuracy  obtained  from  the  variable  step  size 
Simpson's  integration  subroutine  is  much  higher  than  necessary  so  that 
using  a  less  stringent  error  criterion  results  in  less  sample  points  being 
required  and  thus  less  computer  time. 

The  number  of  sample  points  which  the  Simpson's  integration  subroutine 
requires  depends  upon  the  smoothness  of  the  function  being  integrated.  In 
some  cases  very  narrow  high  amplitude  spikes  occur  ir  the  integrand  for  the 
collision  integral  at  the  edges  of  the  integration  interval  due  to  a 
degeneracy  of  the  integrand.  These  spikes  can  be  eliminated  by  not 
'ncluding  the  end  points  in  the  integration  and  only  approaching  the  end 
points  close  enough  to  obtain  the  desired  accuracy.  This  technique  can 


greatly  reduce  the  number  of  sample  points  necessary  in  degenerate  cases 
while  not  noticeably  affecting  the  accuracy. 

Another  factor  which  has  affected  the  smoothness  of  the  integrand 
is  the  use  of  a  tabular  method  to  evaluate  the  first  integrand  of  the 
collision  integral.  This  has  created  a  quantization  of  the  integrand  with 
discontinuous  steps.  It  is  possible  to  get  around  this  problem  by  reversing 
the  order  of  integrations.  A  smoother  integrand  results  and  the  entire 
program  will  take  less  computer  time. 

9.^  Program  for  the  Next  Quarter.  The  collision  integral  computer 
program  will  be  rewritten  to  change  the  order  of  integrations.  The  computer 
time  for  the  entire  program  will  be  checked  in  order  to  keep  it  within 
reasonable  limits  while  still  maintaining  the  desired  accuracy. 
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